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ABSTRACT 


A set of test “modules” for verifying the quality of a CMOS process at the end of 
the wafer fabrication is documented. By electrical testing of specific structures, 
over thirty parameters are collected characterizing interconnects, dielectrics, con- 
tacts, transistors, and inverters. Each test module contains a specification of its 
purpose, the layout of the test structure, the test procedures, the data reduc- 
tion algorithms, and exemplary results obtained from 3-, 2-, or 1.6-micrometer 
CMOS/bulk processes. The document is intended to establish standard process 
qualification procedures for Application Specific Integrated Circuits (ASICs). 
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PREFACE 


The described procedures have evolved through their use in qualifying Applica- 
tion Specific Integrated Circuits (ASICs) fabricated through the USC/ISI MOSIS 
Project. Besides this Process Monitor, test chips for defect, parametric yield, and 
reliability analyses are under development at JPL. 

The layout of the test structures of this Process Monitor is available from JPL 
on magnetic tape in CIF or Calma/GDSII format for selected CMOS/Bulk tech- 
nologies. Presently 1.2-, 1.6-, 2.0-, and 3.0-//m processes are supported using 
MOSIS design rules. The authors invite critical comments and suggestions for 
improvement of this publication. 
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Chapter 1 


Introduction 


This publication describes a. CMOS Process Monitor for wafer level testing at the 
end of the wafer fabrication process. The test structures provide a “snapshot” 
of important process and device parameters and are key to the qualification of 
accompanying Application Specific Integrated Circuits (ASICs). The modular 
design of the structures allows for their placement in available space on ASIC 
wafers, including in saw kerfs. A Process Monitor Test Strip has been fabricated 
in 3.0-, 2.0-, 1.6-, and 1.2-/xm CMOS/Bulk since 1982 through the USC/ISI MO- 
SIS Project and used as a qualification tool for ASICs. The goal of this publication 
is to establish one part of a “workmanship” standard for ASIC qualification. It is 
aimed at design and test engineers with at least an undergraduate education. Pub- 
lications specifying complementary test procedures for defect, parametric yield, 
and reliability analyses are in preparation. 

The present publication contains a description of the procedures for designing, 
testing, and analyzing the following test structures: 


1. Split-Cross-Bridge Resistor 

2. Contact Resistor 

3. MOSFET Capacitor 

4. Field-oxide MOSFET 

5. MOSFET: Gross parameters 

6. MOSFET Quartet: Model parameters 

7. Inverter 

8. Ring Oscillator 1 

9. Timing Sampler 1 
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The process and device parameters that are extracted from each structure are 
listed in Table 1.1. An overview of the test structure layout is given in Figs. 1.1 
and 1.2 by a plot of the CMOS/Bulk Test Strip. All test structures are modular in 
design using 2 x N probe pad arrays. The details of the standard probe pad array 
(No. 1) and a compressed probe pad array (No. 2) are shown in Appendix A. For 
pad identification, a modular counting system is recommended which assigns odd 
numbers (1,3, ..., 2 N — 1) to the lower row of pads and even numbers (2,4, ..., 2N) 
to the upper row counting from left to right. 

The description of each test structure is contained within a test module that 
includes the following elements: 

1. Test module name and structure size 

2. Purpose of the test structure 

List of extracted process/device parameters 

3. Geometrical description/design principles 

List of critical design parameters 

4. Test procedure 

4.1 Circuit diagram 

4.2 Test conditions 

List of test input/output parameters 

4.3 Data reduction algorithms 

List of auxiliary process/device parameters 

5. Test results 

List of extracted parameter values 

6. References and further reading 


In this publication the test module name is given in the chapter heading. The size 
is given in terms of required terminals, i.e., number of probe pads. For structures 
which do not fit in the space between the pads, the additional external space 
requirement is roughly estimated in an equivalent number of pads. 

The purpose of the test module is generally stated by a list of process or device 
parameters which we want to extract. (Intermediate parameters are listed in 
parentheses.) 

The geometrical description is given in terms of iayer/mask names used by 

'These more complex test structures are somewhat beyond the scope of a process monitor. 
Furthermore, their design depends quite on the specific process and its rules. They are, therefore, 
treated here rather briefly, and the interested reader should refer to the cited literature. 
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Table 1.1: Parameters Extracted from the CMOS Process Monitor 


Test Structure 

Parameter 

Parameter Name 

Split-Cross-Bridge 

R„ 

Sheet resistance 

Resistor 

w e ,s e 

Minimum linewidth, spacing 


Pe 

Pitch 


Qp 

Quality factor 

Contact Resistor 

Rr 

Interfacial contact resistance 


P 

Interfacial contact resistivity 

MOSFET Capacitor 

T„ 

Gate oxide thickness 


c„ 

Gate-oxide capacitance/area 


C' g „ 

Gate-overlap capacitance /length 

Field-Oxide MOSFET 

Vt field 

Field oxide threshold voltage 

MOSFET: 

iDSoff 

Channel leakage current 

Gross Parameters 

IpBleak 

Drain diode leakage current 


-f G Bleak 

Gate leakage current 


I D Son 

Transistor on current 


Vl)Bbd 

Drain diode breakdown voltage 

MOSFET Quartet: 

K P 

Intrinsic conduction factor 

Model Parameters 

AT, AW 

Decrease in channel length, width 

(Main parameters, 

8 

Secondary body coefficient 

complete list 

€ 

Velocity saturation coeff. 

in Table 7.1) 

Mb, V 

Gate-field, Body- voltage, Drain- 
field mobility degradation coeff. 


X 

Channel length modulation coeff. 



Body effect, Body-current, M- 
factor body effect coeff. 


M 0 

Intrinsic M-factor 


Vit 

Interface trap voltage 


V F bo 

Flatband voltage less V ft 

Inverter 

^7 inv 

Inverter threshold voltage 


V()L , Vq If 

Inverter output voltages 


G 

Inverter gain 


V n m 

Inverter noise margin 

Ring Oscillator 

- 

Gate delay 

Timing Sampler 

- 

Inverter propagation delay 
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Figure 1.1: CMOS/bulk test strip based on probe pad array No. 1, Part 1 
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Figure 1.2: CMOS/bulk test strip based on probe pad array No. 1, Part 2 


1-5 













DEVICE 

CURRENT 

VOLTAGE 

DC 

AC 

PULSE 

DC 

AC 

PULSE 

MEASURE 

'(A) 

ir(t) 

'p(A) 

v( 


v r(|) 

v p(v) 

FIXED 

SOURCE 

■ 1 $ 

, A 
r Y 

'P^ 

v- 

— 

%(~) 

v p (ri) 

VARIABLE 

SOURCE 

'10 

. A/ 

•r ~ 

Y 

l p <y 

w 

v 7 



VP^ 


NOTE: SUBSCRIPT r IMPLIES AN rms (ROOT MEAN SQUARE 
VALUE) AND SUBSCRIPT p IMPLIES A PEAK VALUE 


Figure 1.3: Electrical Symbols for Current, Voltage, Measure, and Source Devices. 


MOSIS for the P-well process. Conversion to N-well or twin-well is left to the 
reader. MOSIS’ layer definitions and design rules are given in Appendix B. Note 
specifically that ACTIVE means an area under thin oxide, which becomes N + - 
doped (e.g., source/drain of a transistor) except under POLY (e.g., gate of transis- 
tor), where it has the doping of the N-bulk or P-well, and under SELECT, where 
it becomes P + -doped. The geometrical description is accompanied by a list of 
critical design parameters. Some parameters are given in multiples of A, where 2 A 
is the minimum gate length of the process, e.g., A = 0.6 pm for a 1.2-pm process. 

The circuit diagrams contain a symbolic or equivalent circuit representation of 
the device under test connected to symbolized measurement instruments. The 
symbols used in this publication are listed in Figs. 1.3 to 1.4. Where possible, 
the symbols were taken from an IEEE standard [1.1]. Where a symbol could not 
be found, a commonly used symbol is indicated. 

The data reduction algorithms start with the listed test output parameters and the 
auxiliary process/device parameters as input and generate values for the extracted 
process/device parameters listed in the “purpose” section. 

Finally, exemplary test results are given as an orientation to the test engineers in 
judging their own data. 
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CIRCUIT ELEMENT 

SYMBOL 

1. POWER TERMINATION 

T 
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3. WIRE INTERCONNECT 

— 
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— 
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Figure 1.4: Electrical Symbols for Wire Crossing, Interconnect, and Terminations. 
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Figure 1.5: Electrical Symbols for Devices. 
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New York (1987) 
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Chapter 2 

Split-Cross-Bridge Resistor 


2.1. TEST MODULE: The Split-Cross-Bridge Resistor is an 8-terminal structure 
fitting completely within a 2 x 4 pad array. 

2.2. PURPOSE: Split-Cross-Bridge Resistors are used to measure sheet resistance 
as well as width, spacing, and pitch of conducting lines designed to their minimum 
rule. One structure is needed for each conducting layer, i.e., the polysilicon layer, 
the (junction isolated) diffused layers, and the metal layers. The relative deviation 
of the extracted pitch from the designed one serves as a quality factor indicating 
geometry-dependent processing and/or other problems like a magnification error, 
lack of homogeneity in sheet resistance or linewidth, and probe/instrument errors. 
The extracted parameters are listed in Table 2.1. The extracted values of the 
geometrical quantities have been given the index e in order to distinguish them 
from the design values. 

Table 2.1: Parameters extracted from Split-Cross-Bridge Resistor 


Parameter 

Parameter Name 

R s 

Sheet resistance 

(W be 

Bridge linewidth) 

w e 

Minimum design rule linewidth 

S e 

Minimum design rule line spacing 

Pe 

Line Pitch 

Qp 

Quality Factor 


2.3. GEOMETRICAL/DESIGN PRINCIPLES: The Split-Cross-Bridge Resistor 
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[2.1], shown in Figure 2.1, consists of a van-der-Pauw cross [2.2] (left), a bridge 
section (middle), and a split-bridge section (right). The split-bridge is layed out 
so that design linewidth W and spacing S have the minimum dimensions set by 
the design rules. The design line pitch is then P = W ~f S. The design linewidth 
of the (unsplit) bridge must be layed out as Wb — 2W + S. The width of the 
four arms of the cross is usually chosen as Wb, too. Four additional geometrical 


Table 2.2: Critical design parameters for Split-Cross-Bridge Resistor 


Parameter 

Parameter Name 

Value 

W 

Split-bridge linewidth 

Minimum design rule 

s 

Split-bridge spacing 

Minimum design rule 

Wb 

Bridge linewidth 

W + 25 

p 

Line Pitch 

w + s 


rules must be followed when designing the Split-Cross-Bridge Resistor so that the 
results are accurate to better than one percent [2.1]. Rule No. 1 requires that 
the length of each of the four arms of the cross be at least twice the width of the 
arm ( L c > 2 Wb) [2.3]. This allows the use of the van-der-Pauw Equation (2.1) 
for the sheet resistance. Rule No. 2 requires that the voltage taps of the bridges 
(terminals 4 to 7) be of minimum width W and that the distances between the 
taps (Lb and L s ) be large enough so that the as-drawn distance between taps can 
be used directly in the calculations. Rule No. 3 requires that the length of the 
bridge voltage taps be at least twice the tap width. Rule No. 4 requires that the 
separation of a voltage tap from a discontinuity in the line be at least twice the 
linewidth. The critical design parameters are summarized in Table 2.2. 

2.4. TEST PROCEDURE: 

2.4.1. Sheet resistance, R s : 

2.4.1. 1. Circuit diagrams: The circuit diagrams shown in Figures 2.2 and 2.3 
indicate that a current source is applied to two adjacent terminals of the van-der- 
Pauw cross and a voltmeter is applied to the other two terminals. The remaining 
terminals are allowed to float electrically. 

2.4. 1.2. Test conditions: R s is determined from the average of four resistance 
measurements. The first two resistances are determined by measuring V 24 , i.e., 
voltage between terminals 2 and 4, for I 13 , i.e., a current between terminals 1 and 
3, forced in each direction (Figure 2.2). The second two resistances are determined 
analogously by measuring the voltage V :i4 for a current I l2 forced in each direction 
(Figure 2.3). The magnitude of these “van-der-Pauw” currents is chosen the 
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□ METAL 


Figure 2.1: Layout of the Split- Cross-Bridge Resistor. 
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Figure 2.2: Split-Cross-Bridge Resistor sheet resistance, R s , measurement circuit, 
configuration No. 1. (The van-der-Pauw structure has been represented by a 
6- resistor network.) 
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Figure 2.3: Split- Cross-Bridge Resistor van-der-Pauw sheet resistance, R a , mea- 
surement circuit, configuration No. 2. 
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same for all four van-der-Pauw voltage measurements so that only the magnitudes 
of the four measured voltages need to be averaged. This current magnitude, 
called I v , is adjusted so that the magnitude of the measured voltages is between 
1 and 10 mV. This condition allows the voltage to be measured accurately with a 
floating voltmeter having 1 /zV resolution and also prevents spurious current flow 
in junction-isolated layers due to junction forward biasing or junction breakdown. 
A list of test input/output parameters is given in Table 2.3. 

Table 2.3: Test input/output parameters for Split-Cross-Bridge Resistor 


Input parameter 

Value 

Output parameter 

Il3 

±/y 

V»(±I v ) 

L\ 2 


v 3i (±J„) 

1 18 

±h 

Vm(±h), Vs7(±4) 


2.4. 1.3. Data reduction algorithm: R s is derived from the van-der-Pauw Equation 

[ 2 . 2 ]: 


Rg — 


JLYil 

In 2 I v 



( 2 . 1 ) 


In this equation V v — (l^ (+/„)) + \ V 2 .i(—I v )\ + | Km (+/„)! + is the 

average magnitude of the van-der-Pauw voltages with I v — |/ 13 | = 1-^12 1 ? cf. Fig. 2.2 
and Fig. 2.3. 

2.4.2. Effective bridge linewidth, Wt, e ' This is an intermediate parameter needed 
for the calculation of the effective line spacing and pitch, see below. 

2. 4. 2.1. Circuit diagram: See Figure 2.4. 

2.4. 2 . 2 . Test conditions: The bridge voltage is measured between terminals 
4 and 6 while a bridge current J 18 is forced between terminals 1 and 8. The 
magnitude of that bridge current is called R and is set so the measured voltage is 
between 1 and 10 mV. The procedure is repeated with the current reversed and 
the voltage magnitudes averaged. 

2.4. 2. 3. Data reduction algorithm: The effective bridge width Wb e is derived from 

Wi. = ^7 — , (2.2) 

Vb 

where V& = (\V 46 (-\-R)\ + \Vie{—R)\)/^ 1S the average magnitude of the bridge 
voltages and R = |Pi8|> cf. Figure 2.4, and Lb is the designed distance between 
the taps shown in Figure 2.1. 


2,4.3. Effective (split-bridge) linewidth, W e : 
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v 46 



Figure 2.4: Split-Cross-Bridge Resistor bridge linewidth, Wbe , measurement cir- 
cuit. 


2.4.3. 1. Circuit diagram: See Figure 2.5. 

2. 4. 3. 2. Test conditions: The split-bridge voltage V 57 is measured between termi- 
nals 5 and 7 while a split-bridge current J 18 is forced between terminals 1 and 8 
with its current magnitude, I g , set so the measured voltage is between 1 and 10 
mV. The procedure is repeated with the current reversed and the voltage magni- 
tudes averaged. 

2. 4. 3. 3. Data reduction algorithm: The linewidth is 


m = 


2V, 


(2.3) 


In this equation, V s = (| V 57 (+I a )\ + \V 57 (—I s )\)/2 is the average magnitude of the 
split-bridge voltage with I s = |J 18 | , and L a is the distance between the voltage 
taps of the split bridge shown in Figure 2.1. Note the factor 2 in Eq. 2.3 indicating 
that only one half of the external current is flowing through the tapped section of 
the split-bridge. 


2.4.4. Effective line spacing, 5 e , effective pitch, P e , and quality factor Qp : 

2. 4. 4.1. Data reduction algorithm: The effective line spacing is given by, cf. 
Fig. 2.1, 

S e = W be -2W e , (2.4) 
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where W be and W e are given by Eqs. 2.2 and 2.3, respectively. Similarly, the 
effective pitch is 

P e = W be ~W e . (2.5) 

As a quality check, the extracted pitch, P e , is compared to the design pitch, 
P = W + S — W b — W. The quality factor for the Split-Cross Bridge Resistor 
expressed in percentage is 

Qp = x 100 , ( 2 . 6 ) 

where P € is given by Eq. 2.5. 

2.5. TEST RESULTS: Exemplary test results are given in Table 2.4 


Table 2.4: Results from Split-Cross-Bridge Resistors 


Parameter 

Poly 

Metal 

W + -Diffusion 

R, (fl/n) 

16.8 

39.3m 

36.6 

w e ( w ) M 

2.08 (3.00) 

4.36 (4.50) 

5.45 (4,50) 

Se (S) [pm] 

3.92 (3.00) 

4.61 (4.50) 

3.57 (4.50) 

Pe ( P ) [pm] 

6.01 (6.00) 

8.97 (9.00) 

9.02 (9.00) 

Qp [%} 

+0.2 

-0.3 

-0.2 


2.6. REFERENCES: 

2.1. M. G. Buehler and C. W. Hershey, “The Split-Cross-Bridge Resistor for Mea- 
suring the Sheet Resistance, Linewidth, and Line Spacing of Conducting 
Layers,” IEEE Trans. Elec. Dev., ED-33, 1572 (1986) 

2.2. L.J. van der Pauw, “A Method of Measuring Specific Resistivity and Hall 
Effect of Discs of Arbitrary Shape,” Philips Res. Rept., 13, 1 (1958) 

2.3. M. G. Buehler, S. D. Grant, and W. R. Thurber, “Bridge and van der Pauw 
Sheet Resistors for Characterizing the Line Width of Conducting Layers,” 
J. Electrochem. Soc ., 125 , 650 (1978) 
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Chapter 3 


Contact Resistor 


3.1. TEST MODULE: Two structures are offered in this module, the 4-Terminal 
D-Contact Resistor and the 6-Terminal L- Contact Resistor 1 . The 4- and 6- 
terminal structures fit completely within 2x2 and 2x3 pad arrays, respectively. 
A criterion will be given by which it can be judged whether the more economic 
4-terminal structure is sufficient. 

3.2. PURPOSE: To determine the interfacial resistance, Rj, of a contact between 
two conducting layers and calculate the interfacial resistivity, p, by multiplying 
Ri with the area of the contact. In order to establish that p is a true process 
parameter (i.e., independent of area), the resistances of at least two contacts with 
different areas must be measured. The presented methods have been applied to 
metal/poly, metal/ diffusion, and metal/metal contacts. We found that, if properly 
measured, p is usually indeed independent of area, except in undersized contacts. 
However, mean and standard deviation of p, taken over a wafer, can vary widely 
from wafer to wafer. For that reason we chose a method which measures directly 
this critical parameter in contrast to the otherwise often used Transmission-Line 
Tap Resistor, see, e.g., [3.3]. The extracted parameters are listed in Table 3.1. 


*In this publication the labels D and L characterize the shape of the current/voltage taps 
rather than their arrangement, cf. the following sections. The use of these labels in this sense 
is now common in the literature although not in accord with their origin, see, e.g., [3.3]. 
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Table 3.1: Extracted Contact Resistor Parameters. 


Parameter 

Parameter Name 

{Rad 

4-terminal D-contact resistance) 

(Rcl 

6-terminal L-contact resistance) 

Ri 

Interfacial contact resistance 

p 

Interfacial contact resistivity 


3.3. 4-TERMINAL D-CONTACT RESISTOR: 

3.3.1. Geometrical Description/ Design Principles: The layout of the D-Contact 
Resistor is shown in Figure 3.1. The square contact of length d is probed with a 
voltage tap V( and a current tap 7, in each of the two conducting layers (i = 1,2). 
The taps have the same width as the contact except in the vicinity of the contact, 
where they widen to a “flange” allowing for misalignments between the conducting 
layers and the contact window in the isolating dielectric layer. This dogbone shape 
of the taps is the origin of the label D. With the current taps in line, as shown in 
the figure, and flange widths, 8 \ and 8 2 , approaching zero, the horizontal current 
flow is approximately 1-dimensional, and the measured 4-terminal resistance, R 4 d 
approaches the interfacial contact resistance, Rj. 

For finite flange widths, R 4 d is somewhat larger than Rj due to lateral current flow 
in the flanges around the contact. The dogbone design minimizes the flange effect; 
however, the complex geometry has so far prohibited an analytical treatment 
which would allow calculation of Rj from R 4 d- The following criterion for the 
validity of the approximation Rj & R 4D has been derived from simulations [3.1], 
[3.2], [3.3]: For the case where <^i = 8 2 = 8 , the measured contact resistance R 4 £> 
differs from Rj by less than 20 percent, or O.SR 4 o < Rj < R 4 d , provided 


r = 


Rad 

Rsi + Rs2 



(3.1) 


where Rg\ and Rs 2 are the sheet resistances of layers 1 and 2, respectively, and 
7 = 0.65. This relation is plotted in Fig. 3.2. The shaded area is the region in 
which Rf) differs from Rj by less than 20 percent. 

Noting in Eq. 3.1 that R 4 p w Rj — p/d 2 , we obtain the smallest interfacial 
resistivity which can be reasonably measured with this test structure as 


pmin — l{Rsi + Rs 2 ) 8 l ' 2 d°' 8 . ( 3 - 2 ) 

This p m i n should be minimized by designing d and 8 to the minimum dimensions 
which are compatible with the process design rules and the structural requirements 
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Figure 3.1: Layout of a metal-poly 4- Terminal D-Contact Resistor 

of the D-Contact Resistor. We conclude the following set of design rules for the 
4-Terminal D-Contact Resistor: 

1. Design contact width d = d mfn , where d min is equal to the larger of the 
following minimum design rules: contact width, first layer linewidth, and 
second layer linewidth. 

2. Design flange widths Si = 8 2 — S m i n , where S m i n is equal to the larger of the 
minimum design rules for contact overlap in the two layers. 

3. Design the tap widths W\ and W 2 equal to the contact width d. 

If p min is larger than the maximum value tolerable, e.g. in an ASIC, the D- 
Contact Resistor should not be used. For metal-semiconductor contacts use the 
6-Terminal L-Contact Resistor instead. For metal- metal contacts, see [3.5]. If the 
area independence of p is to be established, a second contact resistor with twice 
the minimum contact width is recommended. Table 3.2 summarizes the critical 
design parameters. 


The simulations which lead to formulas 3.1 and 3.2 were made for hypothetical 
perfectly aligned contacts. Fortunately, effects on R A o of misalignments in op- 
posite directions cancel each other to first order, so that the relations hold for 
averages from samples with random misalignments. For individual 4- Terminal D- 
Contact Resistors, the design flange width 6 in those formulas should be replaced 
by an effective flange width of 26 accounting approximately for the worst case of 
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Figure 3.2: Useful parameter zone for 4-Terminal D-Contact Resistors. In the 
shaded area the measured and misalignment corrected 4- Terminal D-Contact Re- 
sistance, R 4 D , equals the interfacial resistance, R/, within 20 percent. For mea- 
sured contact resistances from which misalignment effects have not been removed 
the value of 8 should be replaced by 28. 


Table 3.2: Critical design parameters for 4-Terminal D-Contact Resistor 


Parameter 

Parameter Name 

Value 

d 

Contact length = contact width 

dmrn [ 5 ] 

81 , 82 

Flange widths in layers 1,2 

&min 

W u W 2 

Tap widths in layers 1,2 

d 
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misalignment. 

In general, if the validity criterion, Eq. 3.1, for the 4- Terminal D-Contact Resistor 
is not met, the results from that structure can be stated only in the form Rf < R\o- 
This will usually be sufficient for metal-to-metal contacts. 

3.3.2. Test Procedure: 

3. 3. 2.1. Circuit diagram: See Figure 3.3. 



Figure 3.3: Measurement of 4-terminal D-contact resistance, R 4 jj 

3. 3. 2. 2. Test conditions: R 4 o is determined from the average of two resistance 
measurements. The resistances are determined by measuring the voltage F }2 be- 
tween points Vj and F 2 for current forced in both directions between points Iy 
and / 2 . If the current magnitude, /, is the same in both the forward and reverse 
directions, then the resistance can be determined from the average of the voltage 
measurements. The magnitude of the current is adjusted so that the measured 
voltage is between 1 and 10 mV. In the case of the metal-metal via, the higher 
current density needed to achieve this voltage may cause Joule heating which 
could affect the results. This possibility should be recognized in any test and an 
upper current bound be set for the test. Joule heating can be identified by taking 
resistance measurements at higher current levels and observing deviations in the 
resistance from values taken at lower currents. 

3. 3. 2. 3. Data reduction algorithm: R 4 d is derived from the following equation: 

R<d = y , (3.3) 
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where Vd = (|V) 2 |) = (|l / ] 2 (+/)| + |Fi 2 (~/)|)/2 and I is the current. The interfa- 
cial resistance is estimated as Rj R 4D for T > 0.65. In the case where R cl p- has 
not been misalignment corrected, T must be calculated with an effective flange 
width of 28. For T ^ 0.65 the result should be stated as Tlj ^ • 

With the same restrictions for T the interfacial resistivity is estimated from 

p = R[d? ^ R^jjd?. (3.4) 

Consideration should be given to the error which is introduced in the result by a 
deviation of the true d from the designed contact width. A visual measurement 
of d should be attempted. 

3.4. 6-TERMINAL L-CONTACT RESISTOR: 

3.4.1. Geometrical Description/Design Principles: The approach of the 6- Terminal 
L- Contact Resistor is different from that of the D-Resistor by simplifying the ge- 
ometry to an analytically tractable form at the expense of a somewhat larger 
flange effect. However, the 6-Terminal L-Contact Resistor is restricted to metal/ 
semiconductor contacts. The structure is shown in Figure 3.4. As the resistance in 
the metal layer is small compared to that in the semiconductor layer, Rs 2 < Rsi , 
voltage drops in the metal near the contact can be neglected. Consequently, only 
one voltage tap and one current tap is provided in the metal layer and the exact 
metal geometry is not critical although minimum rules are recommended. In the 
semiconductor layer the width of the voltage or current probing taps has been 
enlarged to W\ — d - f £ 1? giving the taps the shape of straight legs, which we label 
L. Furthermore, there are now four taps in the semiconductor layer, which allow 
current injection from opposite sides. The effect of misalignment of the contact 
window with respect to the semiconductor layer, shown in the figure, can thus 
be eliminated to first order. The design should of course be made with £ u = 
S 12 — £ 13 — 8] 4 = Although the effect of the flange on the measured con- 
tact resistance will be removed analytically, the effect should be minimized by 
design in the first place. As for the L-Contact Resistor holds a relation similar 
to that of Eq. 3.1 (with 7=2.0), the flange effect is minimized by designing d and 
8\ to their minimum rules. Table 3.3 summarizes the critical design parameters 
for the 6-Terminal L-Contact Resistor. Again, for the establishment of an area 
independent p, a second contact resistor with d equaling twice the minimum rule 
is recommended. 

3.4.2. Test Procedure: 

3. 4. 2.1. Circuit diagrams: See Figs. 3.5 and 3.6. 

3. 4. 2. 2. Test conditions: First a current J 16 of magnitude I is applied between 
terminals 1 and 6 in both directions, and the average magnitude of the voltages 
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IMPLANT 

CONTACT 

METAL 


Figure 3.4: Structure of metal-semiconductor 6-Terminal L-Contact Resistor. 
Contact shown misaligned although designed with Sij = S\, j = 1 , ..., 4. Size 
of metal flange exaggerated for clarity of drawing. 


Table 3.3: Critical design parameters for 6-Terminal L-Contact Resistor 


Parameter 

Parameter Name 

Value 

d 

Contact length = contact width 

minimum design rule [,x2] 

Si 

Flange width in layer 1 

minimum design rules 

W t 

Tap width in layer 1 
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Figure 3.5: Measurement of 6-terminal L-contact resistance, i? 6 L> configuration 
No. 1. Semiconductor flange represented by 8-resistor network. Voltages V 2546 
and 1 ^ 5 , 1 6 can be measured sequentially with one voltmeter. 



V 25,36 





V 45,36 


Figure 3.6: Measurement of 6-terminal L-contact resistance, Rql, configuration 
No. 2 
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between terminals 2 and 5 as well as 4 and 5, (|* 25 ,i 6 |) and ( | Vj 5 ,ie|) , respectively 
recorded, cf. Fig. 3.5 and test procedure for the 4-terminal D- Resistor. Then the 
current J 36 = ±7 is applied between terminals 3 and 6, resulting in the voltages 
( I *25,36 1) and (1*45,361), cf. Fig. 3.6. 


3. 4. 2. 3. Data reduction algorithms: First the misalignment effect on the contact 
resistance is eliminated by calculating the average 


— 


( 1 * 25 , 161 ) + ( 1 * 45 , 161 ) + ( 1 * 25,36 1 ) + (1*45, 36 [) 

4/ 


(3.5) 


This is to first order the contact resistance which would be measured if the contact 
had a symmetrical flange of 8 \ — (£n + ^12 *f ^13 + 

In the next step, the parasitic effect of the symmetrized flange is removed from 
R§l with the help of the Thin-Film Model (TFM). The TFM assumes a functional 
dependence as R§i{p, d, Rs\ , 81 ). This function is (iteratively) inverted to calculate 
p or Rj = p/d 2 as described in the appendix of this chapter. 

3.5. TEST RESULTS: Results for some contacts are listed in Table 3.4. The 
auxiliary parameters Rs\ and Rs 2 stem from Split- Cross Bridge Resistor (SXB) 
measurements. For d and 8 the design values were used. The uncertainty in these 
dimensions introduces a relatively large potential error in the extracted interfacial 
parameters. A probably improved estimate of 8 may be obtained by 8 e = 8 — A IF/2, 
where AW = W — W e is taken from a SXB measurement of the layer under 
consideration. The table lists results from 4-Terminal D-Contact Resistors and 
6-Terminal L- Contact Resistors. The results in the table represent averages over 
several chips for which random misalignments were assumed. Of the two D- 
Contact Resistors only the one making contact to N-diffusion has a T > 0.65, i.e., 
large enough for sufficiently error free results. In the case of the contact to N- 
poly, the D-Contact Resistor allows only upper limits of the interfacial parameters 
to be extracted. The results from the 6-Terminal L-Contact Resistor, although 
from a smaller contact and from a different run, show an interfacial resistivity of 
the N-poly contact far below the upper limit extracted from the 4-Terminal D- 
Contact Resistor. It is highly recommended to calculate how uncertainties in the 
dimensions propagate into the results. For the case of the 6- Terminal L-Contact 
Resistor in the table, an error in 8 of 0.2 pm results in an error of about 22% in 
the interfacial parameters. An error in d affects obviously p much more directly 
than Ri. 
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Table 3.4: Results from Contact Resistors 


Run 

M- 

Con- 

tact 

d 

[/zm] 

8 

[^m] 

Ri d 

m 

R<sl 

[fi] 

R . si 

m/a] 

RS2 

[ft/a] 

r 

R, 

[ft] 

P 

[O/zm 2 ] 

770 

MND 

3.0 

1.5 

6.1 


15.5 

35m 

0.90 

» 6.1 

« 55 

770 

MNP 

3.0 

1.5 

2.5 


15.5 

35m 

0.37 

<2.5 

<23 

95F 

MNP 

1.6 

0.8 


4.9 

24.1 

50m 

0.47 

2.3 

5.8 

Contacts: MND = 

metal/N-poly, MND 

= metal/N-diffusion 




3.6. APPENDIX: The Thin-Film-Model (TFM): 


The approach of the TFM is to replace the original contact area with: a) an 
equipotential core of dimension smaller than the contact and embedded in the 
semiconductor layer, with the semiconductor layer being treated as a two dimen- 
sional sheet; and b) a lumped resistor with a value equal to the interfacial contact 
resistance connected to the equipotential core region. After some simplifying 
assumptions, the potential in the two dimensional sheet is solved by conformal 
transformations. The upper layer is treated as an equipotential plane. The model 
was derived in [3.4] for a 4-Terminal L-Contact Resistor; but the derivation shows 
that it is valid for the 6-Terminal L-Contact Resistor without modification. 

In the TFM, the measured contact resistance, is the sum of two parts, the 
intrinsic contact resistance, 77/, and the parasitic resistance, Rji, due to the flange 
in the semiconducting layer. Thus: 

R§l — Ri + Rfu ( 3 - 6 ) 


where 


and 


R 


Rsi 


f i 


R / 


arcosh 


d 2 


1 — t\ -f~ A a 


7 r 


«) 


x 


VFi 


with A., given below. The parameters t and a are determined from 


q = arcosh 


6 — 2 sc 


-f- Vo— T 


arccos 


bx -f b — 2x 
6(1 — x) 


where 


, W l 

6 "5TTX + ’ 


(3.7) 

(3.8) 


(3.9) 


(3.10) 
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(3.11) 

(3.12) 

(3.13) 


A Sj 

' W.-X’ 

Xf ~ X 2A S) 

A = Lt [ coth {£) - t] - 

and Lj is the transmission length given by 



Lt = yJp/Rsu 

(3.14) 

by requiring: 


q(x — t — a) — w(Si + K) 

(3.15) 

and 

C 


q{x = t + a) = *(s l + x,y 

(3.16) 


The inversion of Eq. 3.9 for a given q = q 0 can be achieved iteratively by starting 
with x 0 id = — 1 and calculating successive improved values 


®neui — ® old / \ • (3.17) 

q(Xoid) 

The procedure is terminated when q(x new ) approximates q 0 within 1%. The re- 
sulting values for t and a are then substituted into Eq. 3.8 to determine Rf\. 

Combining Eqs. 3.6 and 3.7: 

Rgl = (3.18) 

As Rf i is a function of p, Eq. 3.18 must be solved iteratively for p. For a measured 
R gl = R 0 , use of the false position rule is recommended with initial values pi c — e, 
where e is less than the minimum expected p , and phi = RqX d 2 and with updates 

Pnew ~ Plo + p \[^c0 ~ ^6 [.(plo)]- (3.19) 

ttGL\Phi) — tiGL\plo) 

If ReiXpnew) < Ro, then set pi a = p new else p hi = p new . This procedure is repeated 
until \Re L (p 

new )--Rbl < rjy where rj is of the order of the measurement error of 
the contact resistance. 
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Chapter 4 


MOSFET Capacitor 


4.1. TEST MODULE: The MOSFET Capacitor is a 4-terminal structure con- 
nected to a 2 x 2 pad array. As the Capacitor does not fit between the pads, a 
total area equivalent to a 2 x 3 standard pad array is required, cf. Fig. 1.1. 

4.2. PURPOSE: To extract the gate-oxide capacitance per unit area, C' OTl and 
the overlap capacitance between gate and source and drain per unit length, C' go 
of MOS transistors. Further data reduction gives the gate-oxide thickness, T ox . 
The two-sided reduction in gate length due to processing, AT, can be extracted 
from two Capacitors with different gate lengths. The list of extracted MOSFET 
Capacitor parameters is given in Table 4.1. The structure can also be used for 
extraction of DC MOS transistor parameters. 


Table 4.1: Extracted MOSFET Capacitor Parameters. 


Parameter 

Parameter Name 

C' ox 

Gate-oxide capacitance/unit area 


Gate-to-source and drain overlap capacitance/unit length 

T 

OX 

Gate-oxide thickness 

A L 

Reduction in gate length 


4.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The MOSFET 
Capacitor, shown in Figure 4.1, is an MOS transistor with a gate in the shape of 
a large ring. The inner and outer radii of that ring are designed as R — L/2 and 
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Figure 4.1: MOSFET Capacitor (closed geometry transistor) 
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R + Zr/2, respectively, where L is the design length of the gate. Metal connectors 
run from source, gate, drain, and body to the four terminal pads. In order to 
reduce parasitic series resistances, the single contacts to source and drain shown 
for simplicity in the figure should in reality be replaced by multiple minimum ge- 
ometry contacts placed around the periphery of the gate in the source and drain 
regions. The total gate area of the MOSFET should be sufficiently large to pro- 
vide at least 1.0 pF total gate-oxide capacitance. The critical design parameters 
are the radius of the gate, R , and the gate length, L. The MOSFET Capacitors 
shown in Figure 1.1 for the 3-/xm CMOS process have the critical dimensions of 
R = 54 fim and X — 6 or 24 pm. The critical design parameters for MOSFET 
Capacitor are shown in Table 4.2 

Table 4.2: Critical design parameters for MOSFET Capacitor 


Parameter 

Parameter Name 

Recommended Value 

R 

Gate radius 

« 40 A 

H i), m 

Gate lengths 

« 5, 15 A 


4.4. TEST PROCEDURE: 

4.4.1. Circuit diagram: In Figure 4.2 the MOSFET Capacitor is modeled by a net- 
work of capacitors, part of which are switched on and off according to different bias 
conditions, as explained below. The source and drain terminals are connected ex- 
ternally and treated from now on as one terminal, labeled SOURCE/DRAIN. The 
capacitances of interest are the gate-oxide capacitance, C ox , and the gate-overlap 
capacitance, C go . The body capacitance, ( 7 *,, the source/drain diode capacitance, 
Cd , the probe pad capacitance, C g b , and the offset capacitance, C 0 //, are parasitic 
capacitances. A DC voltage source, Vgds ? and an AC voltage source, u 9 , are con- 
nected in series to the GATE terminal; an AC ammeter, i m , is connected to the 
SOURCE/DRAIN terminal; and an AC voltmeter, v am , is connected between the 
GATE and SOURCE/DRAIN terminals. 

4.4.2. Test conditions: In the test set-up shown in Figure 4.2 it is assumed that the 
impedance of the AC voltmeter is high enough so that its measurement current 
may be neglected. In addition, the impedance of the AC ammeter is assumed 
to be very low so that the SOURCE/DRAIN terminal is a virtual ground. This 
means that the current through Cd may be neglected. By grounding the BODY 
and measuring the AC signal voltage, v 9m , the effects of ( 7 *,, < 7 ^, and C g b can be 
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Figure 4.2: Measurement of the gate-oxide and gate-overlap capacitances using 
a MOSFET. When the body is inverted ( Vgds = 5 V) switches SI and S2 are 
closed. When the body is depleted (Vgds — 0), only switch SI is closed. When 
the body is accumulated (Vgds = —5 V), both switches are open. 
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neglected. The remaining parasitic capacitance, C Q jf , is independent of the device 
operating conditions and is measured as described below. 

“Measured” capacitances are derived from 

f^meaa /27T f 9 V S rm (^*1) 

where /, is the frequency of the AC source. A typical value for the AC source 
voltage is v a = 100 mV(rms) with a frequency between 100 kHz and 1.0 MHz. 
Commercial capacitance meters can be found which are based on the same prin- 
ciple and provide the necessary circuitry ending in terminals usually labeled HI 
and LO, cf. Fig. 4.2. 

The capacitance between the gate and the source/drain is measured as a function 
of the DC bias gate voltage, Vgds • In the three-point method described below 
one measures the capacitance at three values of Vgds • The first capacitance, C i, 
is measured at Vgds — + V/tp, the second, C 2 , at Vgds — 0, and the third, C3, at 
Vgds = —Vdd, where, usually, V D d = 5 V. Test input/output parameters for the 
MOSFET Capacitor are listed in Table 4.3. 


Table 4.3: Test input/output parameters for MOSFET Capacitor 


Input parameter 

Value 

Output parameter 

f. 

e [0.1, 1] MHz 


V s 

E [10, 100] mV(rms) 

I’m Crneas 

Vgds 

+ Vp £), £s0, — Vdd 

C u Ci,C 3 


4.4.3. Data reduction algorithm: The data reduction algorithms are based on 
the operating characteristics of the MOSFET as shown in Figure 4.3. Here, an 
n-MOSFET is shown operating in strong inversion, depletion, and strong accu- 
mulation. 

The dependence of the measured capacitance on the gate bias is shown in Fig- 
ure 4.4. Also indicated are the three key data points Cj, (7 2 , and C3, which are 
used to extract the desired capacitances. Care must be taken that C 2 is measured 
for the MOSFET biased in the depletion mode. It is suggested that a capacitance 
versus voltage curve be plotted to verify the (weak) depletion at Vgds = 0. (The 
bias voltage for C 2 should be 0.5 to 1.0 V below the onset of inversion, i.e., the 
threshold voltage.) 

When the MOSFET is in strong inversion, the oxide capacitance, C ox , is con- 
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nected to the source and drain terminals through the channel formed under the 
gate as seen in Figure 4.3a; thus, the measured capacitance is 


C\ = C meas (V GD s = +Vdd) — C Q ff + C go + C ox . (4.2) 

When the MOSFET surface is depleted, the oxide capacitance is no longer con- 
nected to the gate but is connected instead to the ground through Cb as seen in 
Figure 4.3b; thus 

C 2 = C meag (V GDS = 0) — C 0 ff 4- C go . (4.3) 

Finally, when the MOSFET is strongly accumulated, the overlap capacitance is de- 
coupled from the measurement terminals by depletion/inversion of the source/drain 
pn junctions as shown in Figure 4.3c, giving 

C 3 = C mea3 (VGDS — ~Vdd) — C 0 ff- (4.4) 


Then C ox and C go are given simply by: 


and 


Co* = C,-C 2 

c go = c 2 - c 3 . 


(4.5) 

(4.6) 


The gate-oxide capacitance/unit area is defined by C' ox = C ox /gate- area. For the 
annular geometry shown in Figure 4.1, the gate area is rc{R + L e / 2) 2 — tt(R — 
L e J 2) 2 = 2irRL (n where L e is the effective gate length and R is the designed mean 
radius of the gate; thus 

C' = Cox . (4.7) 

The effective gate length is L e — L — AL , where L is the as-drawn dimension and 
AL is the two-sided reduction in the gate length due to processing. AL can be 
determined by using two devices with the same radius but different L (//(l) and 
1/(2)), assuming AL and C' ox are the same for both devices. Then from Eq. 4.6 


AL = 


L{l)C OT {2) - L(2)C or (2) 
C ox ( 2) - Co,(l) 


(4.8) 


where C ox ( 1) and C ox ( 2) are the measured oxide capacitances of the first and sec- 
ond device, respectively. As a second choice, the Split- Cross- Bridge Resistor (see 
that chapter) can be used to determine AT; however, somewhat different results 
are expected when the cross section of the gate conductor is not rectangular. As 
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(b) BODY DEPLETED (C 2 = C off + C go ) 



(C) BODY ACCUMULATED (C 3 = C of{ ) 



Figure 4.3: n-MOSFET cross section depicting the three regions of MOSFET 
operation used to extract the gate-oxide and gate-overlap capacitances. 
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ACCUMULATION 


Figure 4.4: Three-point method for measuring the gate-oxide and gate-overlap 
capacitances from an n-MOSFET, where the threshold voltage is approximately 


0.8 V. 
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Table 4.4: MOSFET Capacitor results 


Parameter 

Units 

3/mi CMOS p-well 

1.6^m CMOS n-well 

CL 

fF//mi 2 

0.73 

1.317 


fF///m 

0.074 

0.090 

T 

* OX 

nm 

47.0 

26.2 

R 

(im 

54 

36 

L 

fim 

6,24 

4,8 

A L 

fim 

1 

0.44 

Channel 

type 

n 

n 


a last choice, A L = 0 is assumed. 


The gate-overlap capacitance/unit length is given by C' go = C 90 /gate-perimeter, 
where the gate-perimeter is 27r(i2 + L e j 2) 4- 27t(jR — L e /2) — 4irR, thus: 


n’ _ 

50 4tt R’ 


(4.9) 


Notice that C' go does not depend on AX. 

The gate oxide thickness, T ox , is calculated from C' ox — e 0 /T 0X , giving: 



where e 0 = e(Si02) = 34.515 x 10 6 pF///m. 


(4.10) 


4.5. TEST RESULTS: Typical results for MOSFET Capacitors fabricated by a 
3-/zm CMOS p-well and a 1.6-/mi CMOS n-well process are given in Table 4.4, 
where the 3-/im CMOS p-well results follow from the data shown in Figure 4.4. 
The above n-MOSFET discussion can be extended to a p-MOSFET by simply 
changing the sign of Vqds- 


FURTHER READING: 

4.1. “Product Assurance Technology for Custom LSI/VLSI Electronics,” Report 
for period October 1985 - September 1987, JPL Publication 89-1 (1989) 
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Chapter 5 


Field- Oxide Transistor 


5.1. TEST MODULE: The Field-Oxide Transistor is a 4-terminal structure fitting 
completely within a 2 x 2 pad array. 

5.2. PURPOSE: To assure that the surface between MOSFETs is electrically 
isolated, i.e., that, under normal conditions, no inversion channels are formed 
under field oxide subjected to fields from crossing poly or metal lines. Onset of 
surface conduction is characterized by the field-oxide threshold voltage, Vt field* 
The Field-Oxide Transistor method is applied to regions with n-type and regions 
with p-type doping under the field oxide. The absolute value of Vj field should be 
larger than the supply voltage, i.e., usually 5 V, by a safe margin. 

5.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The layout is 
shown in Figure 5.1. The channel length, L, should be minimum feature size. 
However, the channel width, W, should be much greater than minimum feature 
size, say W = 5 L. The gate metal is designed to overlap the active region by a 
distance L go to allow for misalignment and overetching effects. A recommended 
value for L go is the poly-gate overlap design rule, denoted 0:PAn in Table B2 
(Appendix). Note that this overlap region becomes iV + -doped when the gate is 
metal. In the case of a poly gate, the overlap regions have the channel doping, but 
will normally be inverted long before the surface under the field oxide inverts. This 
is because the gate oxide is much thinner than the field oxide. Design parameters 
for an n-channel Field-Oxide Transistor in 3-/xm CMOS p-well technology are 
listed with test parameters and a typical result in Table 5.1. 

5.4. TEST PROCEDURE: 

5.4.1. Circuit diagram: The circuit diagram shown in Figure 5.2 indicates that the 
test set-up requires a current source for I os and a voltmeter for the measurement 
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Tfrble 5.1; Field-Oxide MOSFET Parameters, 


Extracted Parameter: 


Parameter 

Name 


Vt field 

Field-oxide threshold voltage 

Critical Design Parameters: 


Parameter 

Name 

Value 

L 

Channel length 

Minimum feature 

W 

Channel width 

5 L 

Lgo 

Gate overlap 

Design rule 0:PAn 

Test Input/Output Parameters: 
Input Parameter Value 

Output Parameter 

Vbs 

0 


Vog 

0 


Ids 

1 fiA 

Vbs 

Result: 

Parameter 

Type or Value 


Gate material 

Metal 


Channel 

n-type 


Vt field 

15 V 



of Vt field- The bulk terminal is externally connected to the source terminal and 
the gate terminal to the drain terminal. 

5.4.2. Test conditions: For simplicity, we define the field-oxide threshold voltage, 
Vt fields as the gate voltage at which a current I os — 1 P A flows through the 
channel of a Field-Oxide Transistor with W/L = 5. Actually, Vt field is measured 
with the current source set to I os = 1 p A while Vbs = 0 and Vbs — Vos • The last 
condition ensures that the Transistor is operating in the saturation region. It may 
happen that the gate oxide breaks down before the surface under the field oxide 
becomes inverted. In such a case, the measured voltage represents a lower bound 
for Vp field- If a distinction is needed, the current flowing into the gate should be 
monitored. 
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Chapter 6 


MOSFET Gross Parameters 


6.1. TEST MODULE: The MOSFET Gross Parameters module uses a four- 
terminal enhancement-mode MOSFET fitting completely within a 2 x 2 pad array. 
The test can be applied to any of the four transistors of the MOSFET Quartet 
described in the next chapter. 

6.2. PURPOSE: To evaluate the gross parameters of the MOSFET which are 
listed in Table 6.1. 


Table 6.1: MOSFET Gross Parameters 


Parameter 

Parameter Name 

losoff 

Source- drain leakage current 

iDBleak 

Source and drain diode leakage current 

iGBleak 

Gate leakage current 

Son 

Source- drain on current 

Vdbm 

Source and drain diode breakdown voltage 

VDBGbd 

Source and drain isolation voltage 


6.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The layout of 
a four-terminal n-MOSFET in a p-well is shown in Figure 6.1. The MOSFET 
elements are labeled as follows: 5 for source, D for drain, G for gate, and B for 
body. The source and drain series resistances are designed to have a minimum 
value by minimizing the distance between the contacts and the gate. 

6.4. TEST PROCEDURE: All tests are described for n-MOSFETs. Measure- 
ments of p-MOSFETs are performed by reversing the voltages applied to the 
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Figure 6.1: Layout of four-terminal n-MOSFET in p-well 

devices and measuring the magnitude of the current. 

6.4.1. Source-drain leakage current, Ipsoff : 

Circuit diagram: See Figure 6.2. 

Test conditions: iDSoff is the current that flows in the reverse direction through 
the drain diode with the gate biased so that the channel is turned off. Idso/j 
is measured with the body connected to the source, the gate connected to the 
ground, and Vds = 5 V. 

6.4.2. Source and drain diode leakage current, IoBieak '■ 

Circuit diagram: See Figure 6.3. 

Test conditions: Id Bleak I s the sum of the leakage currents that flow through 
the reverse- biased source and drain diodes to the body. Id Bleak will also include 
leakage currents that flow through the gate oxide that overlaps the source and 
drain junctions. iDBleak is measured with the drain connected to the source, the 
gate connected to the body, and Vdb — +5 V. 

6.4.3. Gate leakage current, iGBieak'- 
Circuit diagram: See Figure 6.4. 

Test conditions: IcBieak is the sum of the currents that flow from the gate to 
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Figure 6.4: Gate leakage current, iGBleak , measurement circuit 
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the body, source, and drain with the gate biased so the channel is on. iGBieak is 
measured with the source and the drain connected to the body and Vqb = 10 V. 

6.4.4. Source-drain on current, Ioson- 
Circuit diagram: See Figure 6.5. 

Test conditions: Id Son is the current which flows between the source and the drain 
with the channel turned on. Ioson is measured with the gate connected to the 
drain, the source connected to the body, and Vos = 5 V. 



5 V 


Figure 6.5: Circuit for measuring source-drain on current, Ioson 


6.4.5. Source and drain diode breakdown voltage, VDBbd, and source and drain 
isolation voltage, VDBGbd 
Circuit diagram: See Figure 6.6. 

Test conditions: is the breakdown voltage measured across both the source 

and drain diodes with 1 microampere forced in the reverse-bias direction through 
the diodes. VoBbd is measured with the source and the gate connected to the drain 
and Idb — 1 ^A. 

The source and drain isolation voltage, VDBGbd , is measured the same way ex- 
cept that the gate is also grounded. VDBGbd is somewhat lower than the source 
and drain diode breakdown voltage because of parallel current flow through the 
gate oxide induced by the fields at the edges of source and drain. If VDBGbd is 
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substantially lower than Vdbm a problem with the integrity of the oxide may be 
suspected. This simple, semi quantitative oxide test is intended to serve as a flag 
only. 



Figure 6.6: Source and drain diode breakdown voltage, and source and 

drain isolation voltage, Vdbg’Mj measurement circuit 


6.5. TEST RESULTS: Tables 6.2 and 6.3 show results obtained from one chip 
with sets of four n- and p-MOSFETs with different gate widths, W , and lengths, 
L, cf. the chapter on MOSFET-Quartet Model Parameters. For ease of orien- 
tation, the data are arranged in the tables as they would appear in a W vs. L 
plane. The transistors were fabricated in 1.6-/xm n-well technology. All transistors 
turn on and off satisfactorily. Ioson is approximately proportional to W/L. The 
ratio Iosoff / Ivson is less than 10~ 6 for the n-channel and 10 -7 for the p-channel 
devices. The leakage and breakdown of the source and drain diodes appear to 
be normal. The source and drain isolation voltages are 3 to 4 V lower than the 
diode breakdown voltages but still comfortably large in comparison to the supply 
voltage. 
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Table 6.2: Gross Parameters extracted from n-MOSFET Quartet 


W = 7.2 fim 

L = 1.6 fim 

L = 4.8 fim 

iDSoff 

64 pA 

74 pA 

I DBleak 

5 pA 

6 pA 

f-G Bleak 

100 pA 

192 pA 

7 DSon 

1.70 mA 

0.70 mA 

VDBbd 

19 V 

18 V 

VoBGbd 

16 V 

15 V 

W = 2.4 fim 

L = 1.6 fim 

L — 4.8 fim 

iDSoff 

1129 pA 

114 pA 

iDBleak 

188 pA 

31 pA 

iGBleak 

104 pA 

69 pA 

I DSon 

0.54 mA 

0.21 mA 

VvBbd 

18 V 

19 V 

VDBGbd 

15 V 

15 V 


Table 6.3: Gross Parameters extracted from p-MOSFET Quartet 


W = 7.2 fim 

L = 1.6 fim 

L — 4.8 fim 

iDSoff 

14 pA 

5 pA 

IDBleak 

22 pA 

20 pA 

IGBleak 

22 pA 

2 pA 

7 DSon 

0.72 mA 

0.23 mA 

V D Bbd 

18 V 

18 V 

VoBGbd 

14 V 

14 V 

W = 2.4 fim 

L = 1.6 fim 

L — 4.8 fim 

iDSoff 

933 pA 

10 pA 

IDBleak 

63 pA 

40 pA 

IGBleak 

39 pA 

107 pA 

I DSon 

0.19 mA 

0.06 mA 

VjJBbd 

18 V 

18 V 

VoBGbd 

14 V 

14 V 
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Chapter 7 


MOSFET Quartet: Model 
Parameters 


7.1. TEST MODULE: The MOSFET Quartet consists of a set of four 4-terminal, 
enhancement-mode MOSFETs, each having a different gate geometry. For the 
CMOS process, separate Quartets are required for n-channel and p-channel tran- 
sistors. Each transistor fits completely within a 2 x 2 pad array. 

7.2. PURPOSE: To extract a set of geometry independent parameters which allow 
description of the behavior of 4-terminal MOSFETs by simple model equations. 
First, geometry dependent parameters are extracted for each of the four MOS- 
FETs operating in either the subthreshold, linear, or saturation regions. From 
these “individual” parameters, “global” geometry independent parameters are 
derived. The approach is designed for easy parameter extraction* from a minimal 
data set. 

The MOSFET model used here is most like SPICE MOS3 [7.1] and is similar 
to CASMOS [7.2] and MOSAID [7.3] models. The MOSFET model equations 
are brought into a form which is linear in the unknown coefficients containing 
the model parameters. Fast linear least-squares methods can then be used for 
the parameter extraction [7.4]. Due to the simplicity of the model, the resulting 
parameters still depend somewhat on the extraction equation and the placement 
of points [7.5]. The following procedures are specified to ensure that the same 
parameter values are extracted independent of the test system or operator. The 
MOSFET model parameters are listed in Table 7.1. 

The MOSFET model pertains to the circuit shown in Fig. 7.1, which shows an 
“intrinsic” 3-terminal transistor with added source and drain series resistances, 
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Jf2, and source and drain to body diodes, making up the “extrinsic” 4- terminal 
transistor. Note that the intrinsic gate and drain voltages, V G > and Vp>, are 
referenced to the intrinsic source S’, whereas the (extrinsic) gate, drain, and body 
voltages, Vg, Vp, and Vg are referenced to the (extrinsic) source S. From Fig. 7.1, 
the relationship between the intrinsic and extrinsic voltages is: 


V G . = V G - I C R (7.1) 

and 

V D , = V D - 2 I C R. (7.2) 

Note further that the (extrinsic) drain current, Ip, is the sum of the channel 
current, which is the intrinsic drain current, and the body current, Ig : 


Id = Ic 4 - Ib • (^* 3 ) 

In the following, equations for the channel current will be given which hold 
separately in the linear and in the saturation region. In the subthreshold region, 
the channel current is zero. For the body current, a solution which holds in 
the subthreshold and in the active region will be used. The separation of the 
cited regions is governed by the relation of gate and drain voltages relative to the 
threshold voltage as indicated in Fig. 7.2. The threshold voltage, Vj-, is, in turn, 
a function of the body voltage: 


Vr = Vbi + — Vg 4- V Tsn , (7.4) 

where 

Hi — 2<£/ + Vfbo + Ht (7.5) 

has been called built-in voltage, with <j>f being the Fermi potential in the body with 
reference to the intrinsic level; Vfbo the flatband voltage due to the workfunction 
difference between the gate material and the bulk silicon and due to oxide traps; 
and the flatband voltage shift Vu caused by interface traps [7.6]. The Fermi 
potential depends on the body dopant density, N , like 

<t>f - V t \nN/Ni, (7.6) 


where iV, is the intrinsic carrier density of silicon and V t = kT/q the thermal 
voltage with k being Boltzmann’s constant, T the absolute temperature, and q 
the electronic charge. Again in Eq. 7.4, 7 is the body effect coefficient given by 


7 


y/2e g qN 


(7.7) 
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:ircuit model with source and drain series 
Intrinsic voltages are denoted with primes 
rimed. 


-3 





Table 7.1: MOSFET Model Parameters. 


Parameter 

Parameter name 


Intrinsic Fermi potential 

0 

Conduction factor 

Kp = fJ>o(sox 

Intrinsic conduction factor 

AW 

Decrease in channel width 

A L 

Decrease in channel length 

w e 

Effective channel width 

L e 

Effective channel length 

V T 

Threshold voltage 

Vto 

Threshold voltage for Vp = 0 

v bl 

Built-in voltage 

1 

Body effect coefficient 

8 

Secondary body effect coefficient 

Do 

Delta coefficient 

O’ 

Intrinsic gate-field mobility degradation coefficient 

{9 

Extrinsic gate-field mobility degradation coefficient) 

®B 

Body-voltage mobility degradation coefficient 

R 

Series resistance 

Rp 

Series resistance per unit width 

v' 

Intrinsic drain-field mobility degradation coefficient 

(l 

Extrinsic drain-field mobility degradation coefficient) 

€ 

Carrier velocity saturation coefficient 

A 

Channel length modulation coefficient 

Ibo 

Active region body current 

Iz 

Body-current body effect coefficient 

M 

Subthreshold slope (M-factor) 

M 0 

Intrinsic M-factor 

'Ifm 

M-factor body effect coefficient 

V it 

Interface- trap voltage 

Vfbo 

Flatband voltage less interface- trap voltage 

c 

Intrinsic Vp-dependent fitting factor 

K 

Extrinsic Vp-dependent fitting factor) 

iBcut 

Cutoff current 

VBcut 

Cutoff voltage 

Po 

Fat FET value of parameter p; p — 8,e,T},9,0B, A 

K\vp 

Width coefficient of parameter p; P — D, E, H,T,TB, L 

Klp 

Length coefficient of parameter p; P = D,E,H,T,TB,L 
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where e g is the dielectric constant of silicon and C r ox is the gate-oxide capacitance 
per unit area. The term V Tan contains the short and narrow channel effects. 

The threshold voltage of a “fat” transistor, i.e., one with dimensions so large that 
V'Tsn can be neglected, is denoted as 


V Tf o = V bt + 1 feh. (7.8) 

The validity of Eq. 7.4 is restricted to not too large body voltages. We found 
reasonable fitting for \V B \ < 2 V. Measurements at |Vjg| = 0, 0.75, and 1.5 V are 
recommended. Note that V B values are negative for normal n-MOSFET opera- 
tion. Although positive for normal p-MOSFET operation, they must be inserted 
negative into Eq. 7.4 as well. 

The drain voltages V D < sa i(V G '), which separate the linear region from the satura- 
tion region (cf. Fig. 7.2) are given by 


VD'sCLt 


2{Vg' ~ Fr)/<*c 

1 + ^1 + 2[(e - 0')/a&c 2 ](F c . - V T ) ’ 


(7.9) 


where a = 1 + b = 1 -|- — Vr), and c — 1 -+■ (d* /2 <i6)(Vg' Vt )* Here & is the 

gate-field mobility degradation coefficient, e is the velocity saturation coefficient, 
and S is the secondary body coefficient given by 


6 = r-^-° + S m (7.10) 

sfih-v B 

where D 0 « 7/4 and where 8 gn contains the short and narrow channel terms. 

In the linear region , i.e., Vo 1 < Vo'sat and Fg?' ^ Vr> the channel current is 

/3[V G '-Vt~(1 + 6)Vd'/2]Vd> (7 n) 

c 1 + 0'(V G > - V T ) 4- CV D ’ + 0 B V B 

Here (3 is the conduction factor and £ 0 is the Vp-dependent fitting factor given by 

c = iy'(l - 2 V Z > 4 - v z 2 ,) + e (2 - V Z ')V Z . + (A - |)(1 - V Z .)V Z , (7.12) 


with V z > = V D ./V D 'sat\ the intrinsic drain-field mobility degradation coefficient 
rf] the channel length modulation factor, A; and the intrinsic gate-field mobility 
degradation coefficient, 6 r . The 9 B accounts for the body voltage dependence of 
the mobility [7.7]. 
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The expression of the channel current for the saturation region , i.e., Vd> > Vo> 8 at 
and Vg> > Vr, is derived from the linear region by substituting Vjyi aat for Vo 1 in 
Eq. 7.9 and reducing Eq. 7.12 to (' = e+ A(1 — Vb'/^D'sai)* resulting in 


P[V G > ~ V T ~ (1 + 6)V D >sat/2]V D >sa t 
1 -f- cVb'aat + t V(V C > - V T ) - A (V D , - ^ D'sat ) 4- QbVb 


(7.13) 


An expression for the body current , /#, which holds in the subthreshold region, 
where V t a < Vj and Ic = 0, as well as in the active region, where Vq > Vt and 
Ic > 0, was derived from Swanson’s subthreshold model [7.8] and the Antognetti 
et al. subthreshold- active region model [7.9] by Blaes et al. [7.10] as 


/Bo(l-exp[-V p /V t ]) 

1 + expHVb- V T )/MV t y 


(7.14) 


This equation contains two parameters which still depend on the body voltage, 
the first being 


I Bo — 


2 ^ 2 "^/ - Vb 


(7.15) 


where 7* has the same physical interpretation as 7 given in Eq. 7.5 but is treated 
as an independent fitting parameter, and the second being 


M 0 + -y m /2^/2<t>f - Vb 
1 - V t /2(2<t> f - V B ) ’ 

where for 7 m holds what has been said about 7 z and 


(7.16) 


M„ = 1 + V h /4>, (7.17) 

is related to the interface trap voltage. 

The validity of Eq. 7.14 is restricted to body currents above some minimum value, 
Incut, due to leakage currents and noise in the transistor and in the instrumen- 
tation. A procedure for determining Iscut will be given in the test procedure 
section. 


The parameters /?, 6 , e, 77, 6 , and A, which describe the behavior of individual 
transistors, are dependent on the effective channel width and length, W e = W — 
A W and L e — L — A L y which differ from the drawn quantities, W and L , by the 
the two-sided shrinkages AW and AT, respectively. 

The conduction factor is given by 

W 

P = Kp y 1 , (7.18) 

"e 
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where K P — PoC' ox is the intrinsic conduction factor with being the zero-field 
channel mobility and C' ox the gate oxide capacitance per unit area. 

The dependence of the other individual parameters p = £, 6 , 77 , 0,0#, A, is modeled 

by 

V = Po + K\yplW e + K lp / L e , (7-19) 

where p t , is the parameter value for a “fat” transistor and K PP and Klw are 
mainly empirical constants. 

7.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The layout and 
design principles of a single MOSFET have been briefly described in the previous 
chapter. In order to extract geometry-independent parameters, four MOSFETs 
[7.4] should be designed with the following dimensions. The minimum geometry 
MOSFET has a channel length-to- width ratio of L m /W m . The other three MOS- 
FETs have ratios of NL m /W m , L m /NW m , and NL m /NW m with a recommended 
N = 3. The layout of this MOSFET Quartet is shown in Fig. 7.3. If not prohib- 
ited by design rules, the contacts should be elongated as shown and not created 
from a series of minimum geometry contacts. This allows a more accurate calcu- 
lation of the source-drain series resistance. The contact length, d, should be kept 
constant as well as the distance between contact and gate, Lcg • Both lengths 
should be minimized. 

7.4. TEST PROCEDURE: 

7.4.1. Circuit diagram: See Fig. 7.4. 

7.4.2. Test Conditions: First the threshold voltages are determined, from which 
the boundaries of the three regions of operation can be estimated, cf. Fig. 7.2. 
Then points (V G , Vo) for the measurement of I D are chosen in each region so that 
simplified equations can be used for easy parameter extraction. 

7.4. 2.1. Threshold Voltage: In the model equations presented in the preceding 
section the threshold voltage of a MOSFET has been implicitly defined as that 
voltage at which the potential at the surface of the silicon bulk equals 2 <f>j. This 
condition is difficult to verify experimentally, especially when <j>j is not known 
a priori. Therefore a procedural approach is used here which is based on the 
maximum slope [7.11] of the \fTo versus Vg curve measured in the saturation 
region with Vo = Vdd* For the extraction of the threshold voltage, the following 
approximations are made: 

1. The body current is ignored by setting I p = 0 so that Io — Ic* 

2. The expression for the channel current Eq. 7.13 is simplified by neglecting the 
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Figure 7.4: Measurement of MOSFET model parameters 
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terms with e, 0, A, and 6 . . 

3. The difference between intrinsic and extrinsic voltages is neglected. 
This leads to the simplified saturated drain current 



(7.20) 


It is assumed now that the preceding relation holds at the maximum slope of 
the measured y/T~o versus Vg curve as illustrated in Fig. 7.5. Then the extracted 
conduction factor is 


A 


= 2 



(7.21) 


and the threshold voltage is 


Vt V Gmax 


27 Dmax 


Pr 


(7.22) 


where VG ma x and I Dmax are the coordinates of the maximum slope point on the 
vTo versus Vq curve. This threshold voltage is taken to be representative for the 
particular geometry and body voltage of the measured transistor. The conduction 
factor, however, has been given the index max to indicate that it depends largely 
on the used procedure. It will be determined more accurately later. 

The point of maximum slope is determined by = d 2 F(Ip)/dV^ = 0, where 
F(Id) = v7d« The following procedure is used to find that point and the value 
of the slope from a minimum set of data points: A set of N = 5 equidistant data 
points, centered at Vg start ~ Vcmax with increments of AVer, is acquired and the 5- 
point smoothed second derivative of F in the center point evaluated, as explained 
below. If that derivative is greater (smaller) than zero, an additional equidistant 
data point at a larger (smaller) Vq is taken, the point at the smallest (largest) 
Vg is dropped, and the second derivative is calculated again. This is repeated 
until the second derivative changes sign. Then the point of maximum slope has 
been bracketed, and no further data are needed. Recommended gate voltages 
are: Vg start, = Vg(Id — 1/^A) and A Vg — 50mV. The method is illustrated in 
Figure 7.6. The first and second derivatives of the jth data set are denoted by 
and Fj 2 \ respectively. An algorithm for calculating F^ and F^ from N > 5 
equidistant points by least-squares fitting to a cubic polynomial is described in 
Appendix C. The results for N = 5, with the ith drain current in the jth data set 
denoted by low and i = —2, —1,0, 1,2, are 


_ + 8f ( f W') - F (Id H 2)\ (7.23) 
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GATE VOLTAGE, V G 


Figure 7.5: Threshold voltage derived from the maximum slope method applied 
to the MOSFET saturation region 
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and 


ff 


#F(Ip) 

dvs 

j^miDj-2) - F(i D j~i) 


(7.24) 


2F(I Dj )-F(I Dj+1 ) + 2F(I Dj+2 )\. 


In the figure, the initial five-point data set, labeled Set 0, has a positive second 
derivative, Fq 2 K Additional data sets, labeled Set 1, Set 2, and Set 3, have to be 
formed by measuring additional Io(V(j) values until the second derivative turns 
negative. The values of Vq and F at that inversion point, Vcmax and F max , 
respectively, are then calculated from two-point linear interpolation using the 
values measured on either side of the zero crossing of the second derivative. In 
Fig. 7.6, these two values of the second derivative are F 2 ; and F 3 Thus, Vcmax 

is calculated from 

171(2) 

V Gm „ = V Gk + (2) k (2) AV g , (7.25) 

F k ~ F k+l 

where the index k refers to the left bracketing data set (Set 2 in the figure) and 
k+1 to the right bracketing one (Set 3 in the figure). The term Vc k denotes the 
center gate bias value of Set k. Likewise, F max is calculated from 

F ma , = F k + (F k+1 - F k ) (VGm 2y G Vnk) , (7.26) 

where F k and F k+ 1 belong to the center values of the bracketing data sets. 

Now, the value of the maximum slope of the F(Id) vs. Va curve is determined 
from the first and second derivatives of the two bracketing data sets using cubic 
interpolation: 

Fill = F^ + F^SVo 

+[ A^ (F *+> - *1") - + 2 Fi 2) )]SV* 

+{-^(Fll\ - Fl } ) + (F™ + F< 2) )]#V3 (7.27) 

where SVq = Vcmax ~~ Vc k ■ Implicit in this equation is the condition: Vc, k < 
V(j max < . 

Now the parameters j3 max and Vj can be calculated from Eqs. 7.21 and 7.22, 
respectively, as 

= 2 (Fill) 2 (7-28) 
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Figure 7.6: Five-point maximum slope method for determining the MOSFET 
threshold voltage 



L4 





and 


V T = V Gmal - (7.29) 

t max 

The procedure is carried out for V B = 0, —0.75, —1.5 V. 

7.4. 2. 2. Cutoff body current, I Bcu t' In the subthreshold region, V G < Vt, the 
channel current is zero in our model; hence, the drain current, Ip, is given by 
the body current I B . For Vb > V t , Eq. 7.14 leads to the following simplified 
subthreshold region current expression: 

log I B = log(7 So - I B ) + VG M yJ" lo S e * ( 7 - 30 ) 

A plot of a measured log Ip vs. V G curve is shown in Fig. 7.7. We define I Bcut 
now as that current at which the first inversion point of that curve occurs when 
the current is lowered. The procedure, which searches for the second derivative 
going through zero, is identical to the one described in the previous section if 
one substitutes log7 D for F, V Gcut for V Grnax , and log I Bc ut for F max , resulting in 
Ibcui = 10 Fmair . The procedure is carried out for V B — 0, —0.75, —1.5 V. 

7.4. 2. 3. Subthreshold Current: The drain current is measured in the subthreshold 
region, i.e., V Gcut < V G < V T . For V B — 0,— 0.75,— 1.5 V and a drain voltage of 
Vp = Vdo ? the gate voltage is chosen as V G — Vj{V B ) — sV t — (k f 6)[Vj(Vb) — 
sV t — V Gcu t{V B )} with k = 1,2, 3,4,5; s — 1 for n-channel; and s = 0 for p-channel 
transistors. This makes a total of 15 points in the subthreshold region, as shown 
in Figure 7.8. 


7. 4. 2. 4. Linear Current: The expression for the channel current, ‘Eq. 7.11, can be 
transformed into a function of the extrinsic voltages, defined in Eqs. 7.1 and 7.2 


as 

0{Vc ~ V T - (1 + S)V d /2]V d 

(7.31) 


1 + 9{V g — Vj) -+• (Vp -f OpVp 

with 

6 = 9' + 2/3 R 

(7.32) 

and 

C = C - PR - 2 8(3R. 

(7.33) 

In the derivation of these relations, terms oc 7£ have been neglected. Furthermore, 
if we choose V D <C V Dgat & ( V G — V T ), we can set Co ~ Vo, c f- Eq. 7.12. It is 
recommended that the drain current be measured for body voltages of V B = 
0,— 0.75,— 1.5 V with the drain voltage be set to Vp — 0.05,0.1,0.2,0.5 V and the 
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Figure 7.7: Vc C ut and iBcut derived from the maximum slope method applied to 
the subthreshold region. In this example Vccut = 0.583 V and Incut = 3.97 nA. 
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GATE VOLTAGE, V G (in unit of V D d) 


Figure 7.8: Subthreshold region extraction points (solid points) for body current 
parameters, cf. Eq. 7.31. Open points are threshold voltages determined from the 
maximum slope method. 
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gate voltage to Vq = Vt(Vb) + (k/4)[V[)p — Lr(Vg)] with k — 2,3,4. This makes 
a total of 36 points in the linear region, as illustrated in Fig. 7.9. 



Figure 7.9: Linear region extraction points (solid points) for drain current pa- 
rameters, cf. Eq. 7.32. Open points are threshold voltages determined from the 
maximum slope method. 


7. 4. 2. 5. Saturation Current: The drain current is measured in the saturation 
region, i.e. Vq > Vj and Vg > Vo sa t ~ {Vg ~ Vr)i c f- Fig. 7.2. It is recommended 
that for body voltages of Vg — 0, — 0.75,— 1.5 V, the drain voltage be set to Vg 
“ 0.SVoo^0.9Vod, I.OVod and the gate voltage to Vg = Vr(Vg) 4- {k/4)[VoD — 
Vt(Vb)\ with k = 1,2, 3, 4. This makes a total of 36 points in the saturation 
region, illustrated in Fig. 7.10. 
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V DRAIN CURRENT, \ZTd~ 




GATE VOLTAGE, V Q (in units of V DD ) 


Figure 7.10: Saturation region extraction points (solid points) for drain current 
parameters, cf. Eq. 7.15. Open points are threshold voltages determined from the 
maximum slope method. 
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7.4.3. Data reduction algorithms: In the following reduction procedures, the 
known parameters are listed after the term “GIVEN”, and the unknown pa- 
rameters are listed after the term “FIND”. Parameters are extracted using a 
least-squares fit [7.12] of the model equations to the data. Thus, the number of 
data points required is at least one more than the number of parameters to be 
extracted. Extracted parameters are enclosed in angle brackets in the equations 
shown below. 

7. 4. 3.1. Subthreshold analysis: GIVEN Id{VgiVb) in the subthreshold region 
as specified by the test procedures, Vj{Vb), V t \ FIND Ib 0 , M ( V# ) . The data 
reduction algorithm is obtained from Eq. 7.14. After linearizing this equation by 
taking logarithms and assuming that Vp >> V ( , the subthreshold drain currrent 
Id ~ Ib is characterized by 

ln/ D = (ln/ Bo (F 8 ))l + (7 .34) 

7. 4.3. 2. Linear analysis: GIVEN Id{Vg, Vd, Vb) as specified by the test proce- 
dures, Vr(Vb); FIND /? , 0, i/, 7. The data reduction algorithm is obtained from 
the channel current expression, Eq. 7.31 with C ~ 7, and Ic = Id — Ib- Further- 
more, the S term can be neglected in the extraction region, resulting in 

Y d = mvG-v T -^)-(e)^(v G -v T )-( v )i D -m^-(e B )^v B , ( 7 . 35 ) 

which is first solved for sets of 12 data points Id{Ygi Vd, Pjj) with Vb — const 
in order to extract 8{Vb) and then solved for body voltage independent /3, 0, 17, 
and 8b simultaneously for all 36 data points of each of the four transistors. The 
8{Vb) of this routine will be replaced later by values extracted from the saturation 
region. 

7. 4. 3. 3. Conduction factor analysis: GIVEN 0(W,L)] FIND j Kp, AW , AT, W e , 
L e . The analysis follows from Eq. 7.18. This equation can be algebraically rear- 
ranged to yield the following system of 4 linear equations: 

f3L = (K P )W - { KpAW ) + (A L)p. (7.36) 

For the following equations the effective gate dimensions are now calculated as 
W e = W - AW and L e = L - A L. 

7.4. 3.4. Gate field mobility and series resistance analysis: GIVEN Kp , 0, L e , W e , 
FIND 0o, Aut? Klt or Rp) R- The source or drain series resistance is contained 
in the extrinsic 0, cf. Eq. 7.32. Assuming that the source and drain series resistors 
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have the same effective length L* between the contact and the channel for all four 
transistors, the series resistance can be written as 

R = Rp/W e = R S L*/W e , (7.37) 

where R p is the source or drain resistance per unit width, and R a the sheet resis- 
tance of the source drain regions. 

The combination of this expression with the (3 expression, Eq. 7.18, and substitu- 
tion into Eq. 7.32 leads to the extraction equation 

9 = (0 o ) + ( K WT )/W e + (2 RpKp + K LT )/L e , (7.38) 

where W e , L e , and Kp come from the preceding conduction factor analysis. For 
L > 1.6/zm we obtained good fits of the current in the linear region by assuming 
Klt ~ 0 and extracting Rp from the last term in Eq. 7.38. For L = 1.2/im, 
however, this procedure did not result in good fits. The results improved sub- 
stantially by calculating Rp from Split-Cross-Bridge and Contact Resistor data 
(see those chapters) and extracting Kpp. The effective source or drain resistor 
length is calculated as L* = Lcg + AL/2 + L t coih.{d/ L t ) with the transfer length 
L f — yjp/ R a , where Lcg is the drawn distance from contact to gate and p and 
d are the resistivity and length of the contact, respectively. For consistency, we 
recommend the latter procedure in general. The series resistance is in any case 
calculated from Eq. 7.37. 

7.4. 3. 5. Body-voltage mobility analysis: GIVEN 0#; FIND Op 0 , K\vtb > Kltb • 
The global parameters of the body voltage mobility degradation coefficient are 
extracted from 

9b = (9 bo) 4- (Kwtb)/ W e -f ( KpTB)/Le • (7.39) 

7.4. 3 . 6 . Threshold analysis: GIVEN V T (V B ), W e , L e , 2 <j> f or C' bx \ FIND V bh 7, 
Vtq, 2 <t> f , K wg , K lg • We extract V bi , 7, K\vgj a nd k lg by fitting Eq. 7.4: 

Vt = (Vto) + {i)y/24>, - V B + (. K wg )/W, c - { K wl )/L„ (7.40) 

The goodness of the fit and even the value of 7 is quite insensitive to the value of 
(f>j, so that a rough idea of the doping will lead with Eq. 7.6 to reasonable results. 
While 3-pm processes used doping densities of N ~ 10 15 /cm 3 leading to 2(f) f % 0.6 
V, the doping density in the channel of 1.6-/tim gate transistors is an order of mag- 
nitude higher leading to 2 <f>f « 0.7 V. In the context of the present publication, 
we recommend using the information about the gate-oxide capacitance per unit 
area, C ox = € OX /T OX , extracted with the MOSFET capacitor (see that chapter), 
together with Eq. 7.7 to calculate a second estimate of N. Because of the loga- 
rithmic dependence of (f>f on N , an iterative procedure including Eqs. 7.40, 7.7, 
and 7.6 converges very rapidly. 
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7.4.3. 7. M-factor analysis: GIVEN M(V B ), 2<£ /? V t \ FIND M 0 , 7 m , V it , V FBo ). 
The data reduction algorithm follows from Eq, 7.16: 


M 


i - i 

2 - V B )\ 


= {Mo) + <7m> 


2 


(7.41) 


The interface trap voltage is calculated from Eq. 7.17 as 


V it = (M 0 - 1 )<f> f , (7.42) 

and the flatband voltage due to workfunction difference and oxide traps is calcu- 
lated from Eq. 7.8 as 

V FBo = 2<f>j + V it — Vbi (7-43) 

for an n-MOSFET, and as 

Vfbo = Vbi — 2 (/>f — Vn (7.44) 


for a p-MOSFET. 


7.4.3.8. /^-factor analysis: GIVEN I Bo (V B ), /?, V t , 2 <f> f \ FIND 7 Z . The data 
reduction algorithm follows from Eq. 7.15: 


-I Bo — {l/z ) 


(3V t 2 


2 - V B 

The values for /? come from the linear region analysis. 


(7.45) 


7.4. 3. 9. Saturation analysis: GIVEN Ip(Vg, V Bl V b ) as specified by the test pro- 
cedure, Ibo( Vb), Vj(V b ), (3 , 0, B b \ FIND 8(V b ), e, A. The data reduction algorithm 
follows from the channel saturation current expression, Eq. 7.13, which holds for 
Vg < Vd + Vt(V b ). For Vq > V F (V B ) -f 3V t we can replace I c by I D — I Bo . The 
values for (3 and 9 come from the linear region analysis. The extraction equation 
is 

Ic[ 1 + Q'{V G ' — V T ) -f- 0 B V B ] — (3(Vg> — Vt — Vb^ at /2)Vb' aa < 

= W^b)) — {^WwsaJc + W(Vd' — Vp'satjlc > (7*46) 

where I c « I d ~Ibo with I Bo = y z (3V t 2 /2yj2<j> f - V B and V D > sat is given by Eq. 7.9. 
As V[)i sat depends also on 8 and e, a few iterations between Eqs. 7.9 and 7.46 are 
necessary, starting with e = 8 — 0 in Eq. 7.9. This system is first solved for the 
three groups of twelve data points with V B = const separately in order to extract 
£(Vb). The parameters e and A, which appear to depend on V B after this step, 
are now consolidated by solving the complete system Eq. 7.46 of 36 equations 
simultaneously with £(V#) as a known variable and only e and A as (body voltage 
independent) fitting parameters. 
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7.4.3.10. Secondary body effect analysis: GIVEN 8(Vg), 20/, W e , L e ; FIND Do, 
Kwdi Kld- The data reduction algorithm follows from Eq. 7.10 as 


6 = (Do) 7 1 + ( Kwd)/W , b + ( K LD )/L e . (7.47) 

- Vb 

7.4.3.11. Drain-field mobility analysis: GIVEN tj , (3, R, £, We, FIND t/ 0 , 
Kwh, Klj{. The global parameters of the drain-field mobility degradation factor 
are obtained from 

V = V + 0* = (vo) + ( K W h) /W e + (K L H)/L e . (7.48) 

This equation has been derived from Eqs. 7.33 and 7.12 with term 80R neglected. 

7.4.3.12. Velocity saturation analysis: GIVEN e, W e , L c \ FIND e 0 , Kwe » 

The data reduction algorithm follows from Eq. 7.19 as 

e - (6 0 ) + (K\y E )/W e + {K LE )/L e . (7.49) 

7.4.3.13. Channel length modulation analysis: GIVEN A, W e , L e \ FIND A 0 , K\\ l, 
K ee . The data reduction algorithm follows from Eq. 7.19 as 

A - (A 0 ) + ( K WL )/W e + (K LL ) / L e . (7.50) 

The the test and extraction procedures are summarized in Table 7.2. 

7.5. RESULTS: The model presented above has been successfully applied to MOS- 
FET Quartets fabricated in 2-, 1.6-, and 1.2-/zm technologies. The extracted 
global, i.e., geometry independent, parameters are quite stable and the goodness 
of fit to the experimental data is only slightly deteriorated at 1.2 (im . We present, 
therefore, a typical 1.2-fim result as the worst case. 

Tables 7.3 and 7.4 list the global parameters for n- and p-channel MOSFETs 
extracted from one Quartet each, located on the same chip. Also, the absolute 
and relative errors are given. Most of the parameters with negative values are 
essentially zero except d E , which must be negative in order to compensate for the 
negative sign of Vg. Large differences in AW and A L between n- and p-channel 
devices have been seen in all fabrication runs. 

Tables 7.5 and 7.6 list the individual n- and p-MOSFET parameters which were 
calculated from the global parameters as a function of (W, L, Vg). The correlation 
between the drain currents calculated from these parameters and those measured 
over the entire linear and saturation regions (all data points in Figs. 7.12 and 7.13 
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Figure 7.11: Maximum slope method for saturation region threshold voltage. In 
this example = 1.094 V, iDmax — 6.30^A, (3 = 99.2/zA/V 2 resulting in 

V T = 0.711 V 
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Table 7.2: MOSFET Parameter Extraction Procedures 


Input variables 

Output variables 

Test procedures 

VGstarl J AVg = 50mV,VR, Vo = V[)D 

Vt(V b ),V Gcu( (V b ) 

Vg = Vt(V b) -»V t - (k/6)[V T (V B ) - sV t 

1 

-Vo OT ,(Vb)];fc = l,2,3,4,5 

) j£>(subthreshold) 

s = (1,0) for (n,p)~channel, Vo = Vdd? Fb 

J 

Vc = V T (V B ) + (fc/4)[Vbp - Vj'(Vg)] 

1 

k = 2,3,4 

) /^(linear) 

Vo = (0.05, 0.1, 0.2, 0.5)V, V B 

J 

V G = Vr(Vfl) + (fe/4)[V OI > - V T (V B )] 

l 

k = 1,2, 3, 4 

) //; (saturation) 

Vo = (0.8, 0.9, 1.0)V, V„ 

J 

Da£a reduction procedures 

/p(subthreshold), Vt(Vb) 

^bo(Vb), M(V b ) 

I D ( linear), Vy(Vs) 


p 

K P ,AW,AL,W c ,L e 

0,Kp,W e ,L e ,[R:,p'] 

00, Kwt , [■K’lt]* Rp, R 

9b 

0bo, Kwtb, Kltb 

Vt(Vb), W e , L e , 2tj>p [C‘ 0 l] 

Vbiy 7 , Pro? [ 2 <^y], iifvvG? 

M(V b ), 2<t>, 

-Pfo? T>n ? > PpSO 

I H0 (V B ),/3Mf 

Tz 

I D ( saturation), I B0 (V B ), V T {V B ),f3, 9 , 0# 

S(V B ),e,X 

S(V b )Mi,W„0 

Do, K wd ,Kld 

p,(3,R,W e ,L e 

Voi K\vh,Klh 

e, IT e , L e 

e l)> KlE 

\W e ,L e 

^05 


Throughout V B = (0, —0.75, — 1.5)V; [] optional; * auxiliary parameters 
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top or bottom plots) is for each {W, L, Vg) expressed in the correlation coefficient 
CC. The average correlation coefficients for the n- and p-channel transistors are 
0.9990 and 0.9994, respectively. 

Figs. 7.12 and 7.13 show, as typical examples, measured data points and calcu- 
lated curves for n- and p-channel transistors the correlation coefficient of which 
happened to coincide with the just mentioned averages. The Id vs. Vp plots 
emphasize the saturation region, the Id/Vd vs. Vs plots the linear region, and 
the log Id vs. Vq the subthreshold region. The data points used in the param- 
eter extraction procedures are filled by a cross. Both figures show good fits in 
all regions; however, the n-channel transistor curves show the larger deviations, 
as expected from the smaller correlation coefficient. No systematic dependence 
of CC on (W, L,Vq) has been seen. Correlation coefficients for transistors with 
I<m > 1.6/xm were consistently greater than 0.9990. 
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Table 7.3: Global Model Parameters for n-Channel Transistors of 1.2 — /xm CMOS 
N-Well Process 


GIVEN: 

Rs(Q/sq) 

= 

80.70 





p( S2 /m 2 ) 

= 

77.59 





Tox(nm) 

= 

21.40 




SUBSTRATE : NSUB ( cm- 3 ) 

= 

2.83E+16 




2<f>f(V) 

= 

0.749 




BETA: 

KP(^A/V 2 )~/joQo 

= 

87.609 

+ 

1.869 ( 

2.13%) 


AW(/vm) 

= 

0.383 

+ 

0.038 < 

10.02%) 


AL(jt/m) 

- 

0.362 

+ 

0.037 ( 

10.30%) 

THRESHOLD: 

: VTf0(V) 

= 

0.779 

+ 

0.059 ( 

7.59%) 


VFB(V) 

= 

0.537 

+ 

0.050 ( 

9.33%) 


T(/V) 

= 

0.606 

+ 

0.036 ( 

5.99%) 


KLG(V*/ym) 

= 

0.133 

+ 

0.021 ( 

16.06%) 


KWG( V*/t/m) 

= 

-0.005 

± 

0.037 ( 

-817.91%) 

DELTA: 

D0(/V) 

= 

0.187 

+ 

0.035 ( 

18.76%) 


KLD(//m) 

= 

0.365 

+ 

0.029 ( 

7.95%) 


KWD( //m) 

= 

0.185 

+ 

0.050 ( 

27.22%) 

THETA: 

00 ( 1/V) 

= 

0.071 

+ 

0.007 ( 

9.79%) 


KLT(/ym/V) 

= 

0.061 

+ 

0.006 ( 

9.70%) 


KWT(/ym/V) 

= 

0.067 

+ 

0.010 ( 

15.41%) 

THETAB: 

0b0( 1/V) 

= 

-0.0230 

+ 

0.0041 ( 

-18.00%) 


KLTB ( /um/V ) 

= 

-0.0061 

+ 

0.0035 ( 

-58.05%) 


KWTB(//m/V) 

= 

-0.0456 

± 

0.0061 ( 

-13.47%) 

ETA: 

H0( 1/V) 

= 

-0.048 

± 

0.001 ( 

-1.91%) 


KLH ( /ym/V ) 

= 

0.039 

+ 

0.001 ( 

1.96%) 


KWH(^m/V) 

= 

0.005 

+ 

0.001 ( 

25.86%) 

EPSILON: 

E0( 1/V) 

- 

0.068 

+ 

0.039 ( 

57.18%) 


KLE(/vm/V) 

= 

0.478 

+ 

0.033 ( 

6.96%) 


KWE(yym/V) 

= 

-0.045 

+ 

0.058 ( 

-130.05%) 

LAMBDA: 

L0( 1/V) 

= 

0.000 

± 

0.002 ( 

816.06%) 


KLL ( /tvm/V ) 

= 

0.038 

+ 

0.002 ( 

6.43%) 

M: 

M0{unitless ) 

= 

1.117 

+ 

0.079 ( 

7.11%) 


Vit(V) 

_ 

0.044 

± 

0.030 ( 

68.06%) 


rm(/V) 

= 

0.810 

± 

0.176 ( 

21.69%) 

IB0: 

Tz(/V) 

= 

3.384 

± 

0.317 ( 

9.38%) 
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Table 7.4: Global Model Parameters for p-Channel Transistors of 1.2 - /xm CMOS 
N-Well Process 


GIVEN: 

Rs ( Q/sq ) 

= 

130.00 





p( S2 fj m 2 ) 

= 

70.72 





Tox(nm) 

= 

21.40 




SUBSTRATE : NSUB ( cm- 3 ) 

= 

1.64E+16 




2<|>f ( V ) 

= 

0.720 




BETA: 

KP(yuA/V 2 ) =//oCo 

= 

29.058 

+ 

0.010 ( 

0.04%) 


AW(/ym) 

= 

0.471 

± 

0.001 ( 

0.12%) 


AL(//m) 

= 

-0.030 

+ 

0.001 ( 

-2.49%) 

THRESHOLD 

: VTf0( V) 

= 

0.967 

+ 

0.025 ( 

2.58%) 


VFB(V) 

= 

-0.253 

+ 

0.021 ( 

-8.42%) 


r(/v) 


0.461 

± 

0.015 { 

3.30%) 


KLG(V*/vm) 

= 

0.148 

+ 

0.015 ( 

10.07%) 


KWG(V*//m) 

SB 

0.001 

+ 

0.015 ( 

1027.45%) 

DELTA: 

D0(/V) 

— 

0.183 

± 

0.021 ( 

11.61%) 


KLD(ywm) 

= 

0.304 

+ 

0.027 ( 

9.02%) 


KWD(/vm) 

= 

0.096 

+ 

0.027 ( 

28 . 53% ) 

THETA: 

00(1A) 

= 

0.135 

+ 

0.008 < 

5.79%) 


KLT(/wm/V) 

= 

0.023 

+ 

0.010 ( 

43.91%) 


KWT(jwmA) 

= 

0.080 

+ 

0.010 { 

12.50%) 

THETAB: 

0b0(lA) 

= 

-0.0760 

+ 

0.0025 ( 

-3.25%) 


KLTB(/umA) 

= 

0.0016 

+ 

0.0032 ( 

204.36%) 


KWTB(/ymA) 

— 

-0.0137 

+ 

0.0032 ( 

-23.04%) 

ETA: 

H0(1A) 

= 

0.020 

+ 

0.001 ( 

3.45%) 


KLH(^mA) 

= 

0.012 

+ 

0.001 ( 

7.85%) 


KWH(//mA) 


-0.019 

+ 

0.001 < 

-4.66% ) 

EPSILON: 

E0(1A) 

= 

0.058 

+ 

0.014 ( 

23.70%) 


KLE(/wmA) 

= 

0.154 

+ 

0.018 ( 

11.74%) 


KWE(/ymA) 

= 

0.008 

+ 

0.018 ( 

223.91%) 

LAMBDA: 

L0(1A) 

= 

0.001 

+ 

0.002 ( 

232.22%) 


KLL(^mA) 

= 

0.091 

+ 

0.003 ( 

3.33%) 

M: 

M0(unitless) 

= 

1.300 

+ 

0.105 ( 

8.09%) 


Vit(V) 

= 

0.108 

+ 

0.038 1 

35.03%) 


Tm(/V) 

= 

0.434 

+ 

0.230 I 

52.88%) 

IB0: 

Tz(/V) 

= 

3.434 

+ 

0.215 ! 

[ 6.26%) 
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Table 7.5: Individual Model Parameters for n-Channel Transistors of 1.2 - fim 
CMOS N-Well Process 


L 

W VB 

(3 

VT 

8 

9' Y\ 

8 

X 

R 

CC 

1.20 

5.40 0.00 

5.24E-04 

0.620 

-0.183 

0.131 0.000 

0.493 

0.045 

42.1 

0.9981 

3.60 

5.40 0.00 

1.36E-04 

0.737 

0.140 

0.077 -0.034 

0.070 

0.012 

42.1 

0.9995 

3.60 

1.80 0.00 

3.83E-05 

0.735 

0.234 

0.043 -0.032 

0.048 

0.012 

149.0 

0.9997 

1.20 

1.80 0.00 

1.48E-04 

0.618 

-0.089 

0.097 0.003 

0.471 

0.045 

149.0 

0.9981 

1.20 

5.40 -0.75 

5.24E-04 

0.837 

-0.246 

0.131 0.000 

0.493 

0.045 

42.1 

0.9979 

3.60 

5.40 -0.75 

1.36E-04 

0.955 

0.077 

0.077 -0.034 

0.070 

0.012 

42.1 

0.9999 

3.60 

1.80 -0.75 

3.83E-05 

0.952 

0.170 

0.043 -0.032 

0.048 

0.012 

149.0 

0.9993 

1.20 

1.80 -0.75 

1.48E-04 

0.835 

-0.152 

0.097 0.003 

0.471 

0.045 

149.0 

0.9986 

1.20 

5.40 -1.50 

5.24E-04 

1.004 

-0.274 

0.131 0.000 

0.493 

0.045 

42.1 

0.9989 

3.60 

5.40 -1.50 

1.36E-04 

1.121 

0.049 

0.077 -0.034 

0.070 

0.012 

42.1 

0.9999 

3.60 

1.80 -1.50 

3.83E-05 

1.119 

0.142 

0.043 -0.032 

0.048 

0.012 

149.0 

0.9995 

1.20 

1.80 -1.50 

1.48E-04 

1.002 

-0.180 

0.097 0.003 

0.471 

0.045 

149.0 

0.9982 

L 

W VB 

IB0 

M 

Temp 

ID0 


0b 



1.20 

5.40 0.00 

6.85E-07 

1.609 

300.0 

1.31E-11 


-0.0378 



3.60 

5.40 0.00 

1.77E-07 

1.609 

300.0 

7.86E-12 


-0.0356 



3.60 

1.80 0.00 

5.01E-08 

1.609 

300.0 

3.54E-11 


-0.0556 



1.20 

1.80 0.00 

1.93E-07 

1.609 

300.0 

9.84E-12. 


-0.0640 



1.20 

5.40 -0.75 

4.84E-07 

1.463 

300.0 

3.36E-11 


-0.0378 



3.60 

5.40 -0.75 

1.25E-07 

1.463 

300.0 

1.43E-11 


-0.0356 



3.60 

1.80 -0.75 

3.54E-08 

1.463 

300.0 

1.52E-11 


-0.0556 



1.20 

1.80 -0.75 

1.37E-07 

1.463 

300.0 

2.05E-11 


-0.0640 



1.20 

5.40 -1.50 

3.95E-07 

1.399 

300.0 

2.60E-11 


-0.0378 



3.60 

5.40 -1.50 

1.02E-07 

1.399 

300.0 

3.35E-11 


-0.0356 



3.60 

1.80 -1.50 

2 . 89E-08 

1.399 

300.0 

3.43E-11 


-0.0556 



1.20 

1.80 -1.50 

1.12E-07 

1.399 

300.0 

1.79E-11 


-0.0640 
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Table 7.6: Individual Model Parameters for p-Channel Transistors of 1.2 — fim 
CMOS N-Well Process 


L 

W 

VB 

8 

VT 

5 

0' n e X 

R 

CC 

1.20 

5.40 

0.00 

1.16E-04 

0.847 

-0.012 

0.138 0.026 0.068 0.075 

81.7 

0.9998 

3.60 

5.40 

0.00 

3.95E-05 

0.927 

0.152 

0.125 0.019 -0.014 0.026 

81.7 

0.9998 

3.60 

1.80 

0.00 

1.06E-05 

0.927 

0.205 

0.081 0.009 -0.010 0.026 

302.9 

0.9998 

1.20 

1.80 

0.00 

3.14E-05 

0.848 

0.041 

0.094 0.015 0.073 0.075 

302.9 

0.9994 

1.20 

5.40 -0.75 

1.16E-04 

1.015 

-0.077 

0.138 0.026 0.068 0.075 

81.7 

0.9993 

3.60 

5.40 -0.75 

3.95E-05 

1.094 

0.087 

0.125 0.019 -0.014 0.026 

81.7 

0.9997 

3.60 

1.80 -0.75 

1.06E-05 

1.095 

0.140 

0.081 0.009 -0.010 0.026 

302.9 

0.9991 

1.20 

1.80 -0.75 

3.14E-05 

1.016 

-0.024 

0.094 0.015 0.073 0.075 

302.9 

0.9998 

1.20 

5.40 -1.50 

1.16E-04 

1.143 

-0.105 

0.138 0.026 0.068 0.075 

81.7 

0.9995 

3.60 

5.40 -1.50 

3.95E-05 

1.222 

0.059 

0.125 0.019 -0.014 0.026 

81.7 

0.9996 

3.60 

1.80 -1.50 

1.06E-05 

1.223 

0.112 

0.081 0.009 -0.010 0.026 

302.9 

0.9989 

1.20 

1.80 -1.50 

3.14E-05 

1.144 

-0.052 

0.094 0.015 0.073 0.075 

302.9 

0.9998 

L 

W 

VB 

IB0 

M 

Temp 

ID0 0b 



1.20 

5.40 

0.00 

1.57E-07 

1.580 

300.0 

3.41E-11 -0.0767 



3.60 

5.40 

0.00 

5.33E-08 

1.580 

300.0 

3.94E-11 -0.0792 



3.60 

1.80 

0.00 

1.44E-08 

1.580 

300.0 

1.66E-12 -0.0850 



1.20 

1.80 

0.00 

4.24E-08 

1.580 

300.0 

4.94E-11 -0.0859 



1.20 

5.40 -0.75 

1.10E-07 

1.495 

300.0 

5.48E-11 -0.0767 



3.60 

5.40 -0.75 

3.73E-08 

1.495 

300.0 

5.42E-11 -0.0792 



3.60 

1.80 -0.75 

1.01E-08 

1.495 

300.0 

3.22E-11 -0.0850 



1.20 

1.80 -0.75 

2.97E-08 

1.495 

300.0 

5.88E-11 -0.0859 



1.20 

5.40 -1.50 

8.97E-08 

1.459 

300.0 

2.72E-11 -0.0767 



3.60 

5.40 -1.50 

3.04E-08 

1.459 

300.0 

4.35E-11 -0.0792 



3.60 

1.80 -1.50 

8.19E-09 

1.459 

300.0 

2.62E-11 -0.0850 



1.20 

1.80 -1.50 

2.42E-08 

1.459 

300.0 

1.11E-11 -0.0859 
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Figure 7.13: Characteristics of p-channel MOSFET with correlation CC — 0.9994 
between model (curves) and data (points). Filled points used for parameter ex- 
traction. Emphasis on linear (top), subthreshold (middle) and saturation region 
(bottom). 
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Chapter 8 


Inverter 


8.1. TEST MODULE: The CMOS/Bulk inverter is a 4-terminal structure fitting 
completely within a 2 x 2 pad array. 

8.2. PURPOSE: To measure the first four parameters listed in Table 8.1 char- 
acterizing a three-piece linear approximation of an inverter transfer curve and to 
calculate a noise margin. As Fig. 8.1 shows, Voh is the inverter output voltage 
Vo for nominal low input voltage, Vj — 0, and Vol is the output voltage for a 
nominal high input voltage, V\ = V DD . The point where output equals input volt- 
age is called the inverter threshold voltage, Vn n , and the magnitude of the slope 
of the transfer curve in that point is called the inverter gain G. The parameters 
extracted from the CMOS/Bulk Inverter are listed in Table 8.1. 

Table 8.1: Parameters extracted from CMOS/Bulk Inverter 


Parameter 

Parameter Name 

Voh 

Output voltage high 

Vol 

Output voltage low 

Vf inv 

Threshold voltage 

G 

Gain 

v 

v nm 

Noise margin 


8.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The layout of a 
p-well CMOS inverter is shown in Figure 8.2. It is recommended that the channel 
length of the n-MOSFET, L n , equal the length of the p-MOSFET, L fn and that 
the channel lengths be designed at the minimum rule. In addition, the channel 
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width of the p-MOSFET, W p , should be 3/2 times the width of the n-MOSFET 
and that W n be designed using the minimum rule. Thus, the conductance of 
the p- and n-MOSFETs will be approximately equal since the electron mobility 
is approximately 3/2 times larger than the hole mobility. The critical design 
parameters for the CMOS/Bulk Inverter are summarized in Table 8.2. 



Figure 8.1: Inverter transfer curve and three-piece linear approximation. 
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Figure 8.2: Layout of a p-well CMOS inverter. 
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Table 8.2: Critical design parameters for CMOS/Bulk Inverter 


Parameter 

Parameter Name 

Value 

L n 

n- channel length 

Minimum 

Lp 

p-channel length 

Ln 

w n 

n- channel width 

Minimum 

w P 

p-channel width 

l'5W n 


8.4. TEST PROCEDURE: 

8.4.1. Output voltage high, Voh > and output voltage low, Vol’> 

8. 4. 1.1. Circuit diagram: See Figure 8.4. A supply voltage of Vpp = 5 V is 
connected to the VDD terminal. A voltage source is connected to the inverter 
input and set to V}, and the inverter output voltage, Vo is measured with a 
voltmeter. 


V DD = 5V 




Vi 


OH 


V OL — 


Figure 8.3: Inverter, Vqjj and V OL , measurement circuit. 


8.4.1. 1. Test conditions: The parameters Voh and V OL are the measured output 
voltages for V} — 0 V and Vj = V DD = 5 V, respectively. 

8.4.2. Inverter threshold voltage, Vt t W 
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V D D = 5V 



Figure 8.4: Inverter threshold voltage, Vrinv, measurement circuit. 

8. 4. 2.1. Circuit diagram: See Figure 8.4. The inverter input is connected exter- 
nally to the inverter output and the output voltage measured. 

8.4. 2. 2. Test conditions: The test is set up so that Vo — Vn nv . 

8.4.3. Inverter gain, G : 

8.4. 3.1. Circuit diagram: See Figure 8.5. Same circuit as for inverter output level 
measurement. 

8.4. 3. 2. Test conditions: The inverter gain G is determined from a two-point 
method in which the inverter output voltage is measured at the inverter input 
voltages Vrinv + AF and Vji nv — AF, where AF = 25 mV. 


Table 8.3: Test input/output parameters for CMOS/Bulk Inverter 

Input Parameter Output parameter 

Output logic levels: 

F/ = 0, Vpp Vo = Vqhi Vol 

Inverter threshold voltage: 

Vi = Vo Vo = Vrinv 

Inverter gain: 

Vr = V T inv± 25 mV Vo = V Ttnv ± Gx 25 mV 
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V DD = 5V 




V 0H 

V 0L 


Figure 8.5: Inverter gain, G, measurement circuit. 


8.4. 3. 3. Data reduction algorithm: The inverter gain is determined from: 


G = 


Vp(Vi = Vrinv - AVQ - VojV! = V Ttnv + AV) 
2AV 


( 8 . 1 ) 


8.4.4. Inverter noise margin V nm : 

8.4.4. 1. Data reduction algorithm: The inverter noise margin is calculated from 
the inverter parameters described above using the so-called maximum square 
method [8.1]. This method implies that the noise margin V nm is defined by that 
DC voltage increment which changes the output state of a long chain of inverters 
or that of a bistable latch when applied in an adverse sense to all inputs simul- 
taneously. The method is illustrated in Figure 8.6 using the three-piece linear 
approximations of transfer curves. Two cases must be distinguished depending 
on the size of Vrinv Assuming Vql w 0 an d Voh Vddi we obtain for the case 
when V T inv > V D d/ 2, 

V nm « (1 - — ){Vdd - Vrinv), ( 8 . 2 ) 

and when Vr tn v ^ W 2, 

Vnn, « (1 - ( 8 . 3 ) 


8-6 



< 


> 







v nm M’-G' 1 ) V T 


V n m~<1-G'V DD - V T . ) 

inv 


inv 


Figure 8.6: Inverter noise margin as derived from the maximum square method 
for two ranges of Vji nv . The sides of the squares give the noise margins, V nm . 
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8.5. TEST RESULTS: The logic levels Vql and Vqh should be very close to 0 
and Vjjd = 5V. Deviations can be caused by bad contacts. For the designed 
transistor width ratio the threshold voltage, Vj inv , should be close to Vdd/2. For 
the inverter gain, G , we have measured values between 10 and 20. 
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Chapter 9 


Ring Oscillator 


9.1. TEST MODULE: Ring oscillator or other AC test structure. 

9.2. PURPOSE: To evaluate inverter gate delay. 

9.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The specific ge- 
ometrical description of this device is dependent upon the specific design rules and 
devices used. In Figure 9.1, the logic design of a typical ring oscillator is shown. 
However, to achieve accurate results, several design rules must be followed: 



Figure 9.1: The ring oscillator, where N is an integer. 
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1. The number of stages in the oscillator must be an odd number (usually a 
prime number) that is large enough to allow each inverter to fully reach the 
high or low state after switching. If the output is not allowed to reach its 
final value before the input begins to change again, the ring oscillator will 
oscillate at abnormally high frequencies. 

2. The output buffer inverter should be designed using minimum geometry tran- 
sistors to add minimum capacitive load to the oscillator. 

3. The output buffer should be provided with separate power and ground so 
the power drawn by the oscillator can be accurately measured. 

4. The circuit should include a NAND gate trigger to allow control of the onset 
of oscillation, preventing the propagation of multiple pulses. 

9.4. TEST PROCEDURE: The actual test procedure for this device is to be left 

to the discretion of the manufacturer. 

9.5. FURTHER READING: 

9.1 M. G. Buehler, “Microelectronic Test Chips for VLSI Electronics,” in VLSI 
Electronics: Micro structure Science , 6, Academic Press, 554-556 (1983) 

9.2 C. Mead and L. Conway, Introduction to VLSI Systems , Addison- Wesley, 
Reading, Massachusetts, 235-236 (1980) 

9.3 N. Sasaki, “Higher Harmonic Generation in CMOS/SOS Ring Oscillators,” 
IEEE Trans . Electron Devices , ED-29, 280-283 (1982) 

9.4 T. W. Houston, “Comments on Higher Harmonic Generation in CMOS/ SOS 
Ring Oscillators,” IEEE Trans. Electron Devices , ED-30, 958-959 (1983) 
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Chapter 10 


Timing Sampler 


10.1. TEST MODULE: CMOS/Bulk timing sampler. 

10.2. PURPOSE: To evaluate inverter pair propagation delay. 

10.3. GEOMETRICAL DESCRIPTION/DESIGN PRINCIPLES: The timing sam- 
pler structure allows the direct measurement of the on-chip inverter propagation 
delay. The circuit schematic diagram of this structure is shown in Figure 10.1. 
This circuit consists of an inverter chain and two Muller C-element latches. The 
latches are connected to the inverter chain at the input and output of an inverter- 
pair (nodes b and b'). The critical design dimensions are the width and length of 
the MOSFETs in the inverters. 

10.4. TEST PROCEDURE: Delays are measured for rising and falling step inputs. 
The latch connected to node b will trip when the step input signal, IN, reaches 
node b. The time of arrival of the signal at node b is determined by adjusting 
the step enable signal, EN, to coincide with the time of arrival of the IN signal 
at node b. The inverter-pair delay is determined by subtracting the signal arrival 
times measured at node b from that measured at node b'. 

10.5. FURTHER READING: 

10.1 B. R. Blaes, M. G. Buehler, and Y-S Lin, “Propagation Delay Measurements 
from a Timing Sampler Intended for use in Space”, IEEE Trans. Nucl. Sci ., 
NS-34, 1470-1473 (1987) 

10.2 B.R. Blaes and M.G. Buehler, “Inverter Propagation Delay Measurements 
Using Timing Sampler Circuits”, Proc. IEEE 1989 Int. Conf. Microelec- 
tronic Test Structures , 2, 227 (1989) 
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Appendix A 

Currently Used 2 x N Probe Pad Arrays 


Two probe-pad configurations used in the semiconductor industry are shown in 
this appendix. Probe Pad Arrays No. 1 and No. 2 are shown in Figures A1 and 
A2, respectively. Array 1 is bondable, array 2 is not; however, there is a significant 
space savings when this layout is used. The layout of the test strip using Probe 
Pad Array No. 1 is illustrated in Figure 1.1. The stacking box (indicated by the 
dashed lines) is used when butting multiple probe-pad-array blocks together to 
assure uniform pad spacing and test strip probe-ability. The bounding box, which 
defines the boundary of usable and unusable space, is indicated by the dotted 
lines. Features cannot be located outside the bounding box unless they extend 
beyond the end of the probe-pad contact area such as in the case of the area 
MOSFET and ring oscillator test structures illustrated in Figure 1.1. 
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STACKING BOX 


Figure A2: Probe-pad array No. 2. 
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Appendix B 

Geometrical Design Rules 


The notation for the rule set used to design the CMOS Process Monitor is listed 
in Table Bl. The design rules are listed in Table B2. It is a modification of re- 
vision six of the MOSIS CMOS scalable rule set which is taken from the MOSIS 
CMOS-Bulk design rules, dated February, 1988. This is a double-metal single- 
poly enhancement MOSFET process which does not use split butting or buried 
contacts. The complete design rule set is available from MOSIS/ISI. The design 
rule set given here is for a 2/i.m process. All dimensions used in Table B2 are 
in micrometers. Table B3 is provided to define the process and CIF or CALMA 
GDS-2 layer names. 


A personalized CMOS Process Monitor can be obtained by entering the data 
requested in Tables B2 and B3. A column has been left for the fabricator to 
insert the dimensions of his rules in Tables B2 and B3. The submission of your 
rule set dimensions and layer names to JPL will allow preparation of a personalized 
tape containing the CMOS Process Monitor. When providing your rule set, please 
give any additional information that is relevant to making a successful transition 
of the CMOS Process Monitor to you. That is, if there are additional rules that 
should be known, please list them. 
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Table Bl: Geometrical Design Rule Notation. 



Layer 

Layer 

Dimensions 

Names 

Modifiers 

S = Space 

W = Well 

n = n-type 

0 = Overlap 

A = Active 

p = p-type 

W = Width 

P = Polysilicon 


L = Length 

I = Select 
C = Contact Cuts 
M - Metal # 1 
M2 - Metal # 2 
V - Via 
G = Glass 
XT = MOSFET 
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Table B2: Geome t rical Design Rule Specification. 


Layer 

Symbols 

JPL 

Dimensions 

(/xm) 

Fabricator 

Dimensions 

(fim) 

Symbol Description 

L:nXT 

= 2.0 


n-MOSFET minimum drawn chan- 
nel length 

W:nXT 

= 3.0 


n-MOSFET minimum drawn chan- 
nel width 

L:pXT 

= 2.0 


p-MOSFET minimum drawn chan- 
nel length 

W:pXT 

= 3.0 


p-MOSFET minimum drawn chan- 
nel width 

S:WW 

= 6.0 


Spacing Well to Well (same poten- 
tial) 

W:W 

= 10.0 


Linewidth of Well 

S:AA 

= 3.0 


Spacing Active to Active 

S: AW 

= 5.0 


Spacing Active to Well Boundry 

W:A 

= 3.0 


Linewidth of Active 

W:P 

= 2.0 


Linewidth of Poly 

S:PP 

= 2.0 


Spacing Poly to Poly 

S:PA 

= 1.0 


Spacing Field Poly to Active 

0:PA 

= 2.0 


Overlap Poly over Active 

0:IA 

= 2.0 


Overlap Select beyond Active 

W:I 

= 2.0 


Linewidth of Select 

S:II 

= 2.0 


Spacing Select to Select 

S:IA 

= 2.0 


Spacing Select to Active 

OTA 

= 2.0 


Overlap Select to Active 

0:IXT 

= 2.0 


Overlap Select to MOSFET Gate 

W:C 

= 2.0 


Width of Contact 

S:CXT 

= 2.0 


Spacing Contact to MOSFET Gate 

0:M1C 

= 1.0 


Overlap of Metal 1 to Contact 

0:PC 

= 1.0 


Overlap of Poly to Contact 


continued 
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Table B2: Geometrical Design Rule Specification - continued. 


Layer 

Symbols 

JPL 

Dimensions 

(fim) 

Fabricator 

Dimensions 

(fim) 

Symbol Description 

0:AC 

= 1.0 


Overlap of Active to Contact 

S:CC 

= 2.0 


Spacing Contact to Contact 

S:CA 

- 3.0 


Spacing Contact to p-Active 

W:V 

= 2.0 


Width of Via 

0:M1V 

= 1.0 


Overlap 1st Metal beyond Via 

0:M2V 

= 1.0 


Overlap 2nd Metal beyond Via 

W:M 

= 3.0 


Width of 1st Metal 

S:MM 

= 3.0 


Spacing of Metal to Metal 

W:M2 

= 3.0 


Width of 2nd Metal 

S:M2M2 

= 4.0 


Spacing of 2nd Metal to 2nd Metal 


* = p-Well or twin-well process only 
** = n- Well or twin- well process only 


B-4 



Table B3: Process Definition and Layer Naming Conventions (MOSIS) 

Process (check one): 

p-Well n-Well twin Well 

Contacts currently used in the Process Monitor: 

The Process Monitor uses only three types of contacts. The Cuts 
layers specify contacts to Active and Poly from First Level Metal. 
Contacts from First to Second-Level Metal are called out by the Via 
layer. Direct contact from the Second-Level Metal to Active or Poly 
is assumed to be unsupported unless otherwise specified. 


Layer Name 

CIF 

CALMA Fabricator 

P-Well 

CWP 

53 

N-Well 

CWN 

42 

Active 

CAA 

43 

Poly 

CPG 

46 

P- Select (jP + I mplant) 

CSP 

44 

N- Select (N + Implant) 

CSN 

45 

Contact Cut: Metal to Active 

CCA 

48 

Contact Cut: Metal to Poly 

CCP 

47 

Metal 1 

CMF 

49 

Metal 2 

CMS 

51 

Via 

CVA 

50 

Glass/Passivation 

COG 

52 

Ot.Vipr ( 

) 


Ot.hpr ( 

) 


Other ( 

/ 

-) 



Note that each of these layers translates directly to the corresponding 
mask layer. No Boolean or sizing functions are used. If your software 
requires layers to be defined differently, please explain on a separate 
page. 
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Appendix C 

Smoothed First and Second Derivatives 


If we have a set of N = 2m + 1 data points (xi,yi) y i = — rra, ...,rrc, with equidistant 
Xi , then the first and second derivatives about the center value, x 0 , denoted F W 
and F^ 2 \ respectively, can be calculated by 


fi l ) — 




CjimVi 

MtmAx 


and 


FW = 


d 2 y 
dx 2 


= E 


CrtmUt 

M 2m Ax*’ 


where Ax is the increment between successive Xi values, C lim and C l2 m are the first 
and second derivative smoothing coefficients, and Mi m and M 2rn are normalization 
factors. The smoothing coefficients and normalization factors for the first and 
second derivatives are listed in Tables Cl and C2, respectively. These coefficients 
were derived from least-squares fitting of a cubic polynomial to the data sets [Cl], 
[ 02 ). 

REFERENCES 


Cl A. Savitzky and H. J. E. Golay, “Smoothing and Differentiation of Data by 
Simplified Least Squares Procedures,” Analytical Chemistry , 36, 1627 (1964) 

C2 J. Steinier, Y. Termonia, and J. Deltour, “Comments on Smoothing and Dif- 
ferentiation of Data by Simplified Least Square Procedure,” Analytical Chemistry , 
44, 1906 (1972) 
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Table Cl: Smoothing coefficients for calculating the first derivative at the 
central point o f a set of data least-squares-fitted to a cub ic polynomial. 


N 

13 

11 

9 

7 

5 

m 

6 

5 

4 

3 

2 

M\ m 

24024 

5148 

1188 

252 

12 

C_61 m 

1133 





C_5i m 

-660 

300 




C 41m 

-1578 

-294 

86 



C-31m 

-1796 

-532 

-142 

22 


C-21m 

-1489 

-503 

-193 

-67 

1 

C-Um 

-832 

-296 

-126 

-58 

-8 

C —01m 

0 

0 

0 

0 

0 

C\t m 

832 

296 

126 

58 

8 

@21 m 

1489 

503 

193 

67 

-1 

@31 m 

1796 

532 

142 

-22 


C 41 m 

1578 

294 

-86 



C 51 rn 

660 

-300 




m 

-1133 






Table C2: Smoothing coefficients for calculating the second derivative at the at 
the central point of a set of data least-squares-fitted to a cubic polynomial. 


N 

13 

11 

9 

7 

5 

m 

6 

5 

4 

3 

2 

M2 m 

1001 

429 

462 

42 

7 

C— 62m 

22 





C*-52 m 

11 

15 




C-42m 

2 

6 

28 



C-32 m 

-5 

-1 

7 

5 


C-22 m 

-10 

-6 

-8 

0 

2 

C*— 12m 

-13 

-9 

-17 

-3 

-1 

Co2m 

-14 

-10 

-20 

-4 

-2 

Ci2m 

-13 

-9 

-17 

-3 

-1 

C22m 

-10 

-6 

-8 

0 

2 

@32 m 

-5 

-1 

7 

5 


C , 42m 

2 

6 

28 



C 52 m 

11 

15 




CVlm 

22 
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